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We  have  derived  a linear  kinetic  theory  of  parametric 
instabilities  in  an  infinite,  homogeneous,  fully  ionized 
plasma,  with  no  magnetic  field.  Our  derivation  relaxes 
three  assumptions  usually  made  in  the  linear  theory  of  cer- 
tain parametric  instabilities: 

1) .  We  do  not  assume  that  k , the  wavenumber  of  the 
pump  (the  driving  field  of  the  instability) ,is  much  smaller 
than  any  other  inverse  distance  in  the  plasma; 

2)  . We  do  not  restrict  ourselves  growth  of  the  in- 
stability in  one  dimension;  and, 

3) .  We  consider  a standing-wave  pump,  in  addition  to 
the  usual  travelling-wave  pump. 

Using  this  generalized  formalism,  we  have  studied  the 
important  problem  of  the  parametric  instabilities,  in  three 
dimensions,  pumped  by  a monochromatic,  longitudinal  pump 
of  arbitrary  wavevector.  By  solving  a dispersion  relation 
derived  for  these  instabilities,  numerically  and  analyti- 
cally, we  have  found  a number  of  new,  three-dimensional 
effects : 

1).  Only  under  quite  restrictive  conditions  on 
does  the  oscillating  two  stream  (OTS)  instability  have  its 
fastest  growing  daughter  waves  with  wavenumber",  that  are 
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parallel  to  K^-  when  kQ  k,  these  dauqhter  waves  can  grow 
at  as  much  of  an  angle  to  kQ  as  20°. 

2) .  The  Langmuir-pumped  electron-ion  decay  (EID)  in- 
stability, for  arbitrary  k^,  has  regimes  in  which  it  is  a 
forward  scattering  instability  (as  well  as  the  usual  back- 
scattering  regimes  in  the  case  that  kQ  > ^ 20C/3). 

3) .  A new  instability,  which  is  only  seen  in  three 
dimensions,  is  possible.  This  instability,  the  stimulated 
modulational  (SM)  instability,  can  produce  daughter  waves 

at  an  angle  to  kQ,  with  wavenumbers  slightly  larger  than  kQ. 

We  have  applied  these  results  to  three  physical 
situations  modelled  on  am  electron  stream  quasilinearly 
generating  a Langmuir  spectrum  of  sufficient  intensity  to 
act  as  a pump: 

1) .  Type  III  solar  radio  bursts.  We  find  that  three 
dimensional  effects  are  critical  in  understanding  the  re- 
laxation of  Langmuir  waves  at  distances  less  than  450  R . 

© 

2) .  Auroral  streams.  Our  results  show  a one-dimen- 
sional treatment  is  sufficient  in  this  case,  although 
three  dimensional  effects  are  very  likely  to  be  important 
in  determining  the  saturated  state  of  the  instability. 

3) .  Laboratory  electron  beams.  We  conclude  here 
that  parametric  effects  dominate  the  wave  dynamics  for  a 
"slow"  beam,  but  that  this  is  not  the  case  for  faster 
beams . 

In  none  of  these  cases  do  we  find  parametric  processes  that 
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greatly  affect  the  quasilinear  decay  of  the  beam,  contrary 
to  the  findings  of  several  authors  who  have  examined 
similar  applications  of  parametric  instabilities. 

We  have  also  studied  the  effects  of  a standing-wave 
pump  of  arbitrary  k^,  for  two  parametric  instabilities,  the 
2upC  and  CID  instabilities.  The  most  important  result  of  a 
standing-wave  pump  for  these  instabilities  is  a reduction 
of  the  effective  power  of  the  pump.  This  is  a result  of 
different  modes  of  the  plasma  being  excited  by  each  of 
the  two  components  of  the  standing-wave  pump;  these  ex- 
cited modes  are  coupled  more  strongly  (raising  the  effective 

pump  power)  as  the  (k  -dependent)  frequency  separation  of 

o 

the  modes  decreases  relative  to  the  growth  rate  of  the  in- 
stability. As  kQ  goes  to  zero,  and  a standing-wave  becomes 
indistinguishable  from  a travelling-wave,  we  recover  the 
travelling-wave  results.  When  we  apply  our  results  for  the 
EID  instability  to  ionospheric  modification  data,  we  find 
the  decrease  in  the  effective  pump  power  resulting  from 
a standing-wave  pump,  lowers  the  theoretical  predictions 
for  the  Langmuir  turbulence  produced  in  these  experiments. 
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CHAPTER  I 


INTRODUCTION 

The  coupling  of  normal  modes  of  a linear  system  by 

a time-dependent  interaction  is  a phenomenon  first  noted  by 

Faraday*  in  1831  and  studied  theoretically  by  Lord  F.ay- 
• 2 

leigh  in  1883.  Now  called  parametric  excitation,  the 
study  of  these  processes  provides  a general  and  powerful 
framework  for  studvinq  linear  responses  of  a svstem.  In 
this  dissertation,  the  theory  of  the  parametric  excitations 
of  an  infinite,  homoqeneous , fullv-ionized  plasma  is  ex- 
tended to  include  two  effects  dependent  on  the  wave- 
vector,  1^,  of  the  driving  field  o'  the  instability: 

1) .  The  effects  of  this  drivinq  field  (the  "pumo") 
being  a standing-wave  field  rather  than  the  usual  travellino- 
wave  field. 

2)  . Three  dimensional  effects  of  a longitudinally 
polarized  pump  with  arbitrary  wavevector. 

We  have  found  that  in  a larqe  number  of  physical  situations 
the  spatial  dependence  of  the  parametric  coupling,  depend- 
ent on  k^,  is  imnortant  in  determining  both  qualitative  and 
quantitative  features  of  this  miasma  interaction. 

The  first  systematic  observations  of  parametric  ex- 
citation were  done  by  Michael  Faraday,  when  he  noted  that 


he  could  maintain  water  waves  at  half  the  frequency  of  the 
vibrating  support  which  generated  them.  In  1859,  F.  Melde^ 
published  results  of  an  exhaustive  study  he  had  carried  out 
of  the  excitation  of  transverse  modes  of  a vibrating  wire 
which  could  be  amplified  by  a longitudinal  vibration  of  the 
wire  at  twice  the  transverse  frequency.  Lord  Rayleigh  made 
several  theoretical  studies  of  this  problem,  and  his  re- 
sults exhibited  the  two  characteristic  features  of 
parametric  excitation:  the  existence  of  approximate  fre- 

quency matching  (in  the  examples  given  above,  the  excited 
disturbances  occurred  at  half  the  frequency  of  the  driver, 
or  "pump"),  and  a threshold  power  for  the  pump,  above  which 
the  normal  modes  of  the  system  begin  to  be  amplified. 

A much  more  commonplace  example  of  this  phenomenon  is 
the  energy  transfer  that  occurs  in  a child's  pumping  a 
swing.  The  time-varying  element  in  this  otherwise  simple 
pendulum  is  the  raising  and  lowering  of  the  effective 
length  of  the  pendulum;  the  child  raises  his  center  of 
gravity  as  the  swing  goes  up  and  lowers  it  as  the  swing 
comes  down.  The  frequency  of  this  variation  is  twice  that 
of  the  natural  frequency  of  the  swing.  This  "pumping" 
couples  together  the  two,  here  degenerate,  normal  modes  of 
the  pendulum,  with  the  frequency  matching  condition: 

= I I +•  I / = { dV  [ + \~0j\  - 2 Oj  ■ 

And,  if  this  pumping  is  of  sufficient  amplitude,  then  both 
normal  modes  become  unstable  and  energy  is  transferred  to 
these  modes. 
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In  general,  parametric  excitation  occurs  when  several 

otherwise  independent  resonances  of  a system  are  coupled 
together  by  some  time-varying  parameter  of  the  system.  If 
amplification  of  these  normal  modes  is  to  occur,  the  time 
variation  of  this  parameter  is  not  arbitrary,  but  rather  is 
fixed  by  a condition,  dependent  on  the  frequencies  of  the 
normal  modes  which  are  being  amplified.  When  there  is  dis- 
sipation in  the  uncoupled  modes,  then  a "threshold"  power 
of  the  time-varying  parameter  exists,  above  which  energy  is 
transferred  to  the  now-coupled  normal  modes.  Notice  that 
we  have  made  the  implicit  assumption  that  the  energy  trans- 
fer is  small  enough  that  the  pump  remains  essentially  un- 
affected by  the  excitation  of  modes  in  the  system;  this  is  the 
origin  of  the  term  "parametric,"  the  time-varying  parameter 
being  treated  as  having  constant  amplitude  over  many 
periods , and  entering  the  calculations  as  a parameter  in 
the  growth  rates  of  the  excited  modes.  Current  application 
of  this  principle  is  being  made  in  many  fields,  including 
circuit  design*,  non-linear  optics^,  and  several  areas  of 
plasma  physics. 

In  the  case  of  a plasma,  parametric  couplings  occur 
in  amazing  variety,  as  a result  of  the  multitude  of  waves 
possible  in  a plasma  and  the  richness  of  the  interactions 
among  them.  If  we  restrict  ourselves  to  an  infinite,  homo- 
geneous, unmagnetized,  fully-ionized  plasma,  there  are  only 
three  linear  resonances  possible,  each  with  a well-known 
dispersion  relation: 
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1).  the  electromagnetic  wave,  e-k  * 0 (see  page  12  for 


our  notation)  : 

~ [l  + (c'^AV koi )]  1 k^t.  = 4/rnc2/r v\t- 

6spt 

A - 

2) .  the  Langmuir  wave,  e‘k  = 1, 

SI  - (/♦  sV/ki)54 , kL  = titne'A, ; 

a.  S'  . 

3) .  the  ion-acoustic  wave,  e-k  = 1, 

tjA  f /.  /&&  , &c  , 

‘*1'  (m.An.)*  k/fe*  , ■ 

V. 

Each  of  these  waves  has  a well-understood  dissipation 
which  we  list  here  as  the  imaginary  part  of  the  above  fre- 
quencies : 

1) .  the  electromagnetic  wave,** 

&L  ~ ku.  in(^£*W/[6yinrrK,^c]j 

2) .  the  Langmuir  wave,** 

3) .  the  ion-acoustic  wave,7 

ink*.,  / 0AI  ,6<*6,: 

<*U  | (rfe/tn, ?*{?/*  /*■  , 0 

In  a quiescent  plasma,  these  three  waves  are  independ- 
ent; however,  an  intense  high  frequency  field  can  paramet- 
rically couple  together  various  combinations  of  these 
modes.8  Physically,  this  high-frequency  field  acts  as  the 
pump,  by  setting  up  currents  in  the  plasma.  These  currents 
change  the  susceptibility  of  the  plasma,  nonlinearly,  in 
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such  a way  that  the  equations  describing  the  resultant 
fields  in  the  plasma  have  source  terms  which  depend  on  a 
power  series  in  the  pump  power.  Thus,  for  example,  the 
equation  describing  the  amplitude  of  a low-frequency  ion- 
acoustic  wave  contains  a term  due  to  the  beating  of  a high- 
frequency  pump  and  a high-frequency  daughter  field.  Given 
that  the  appropriate  frequency -matching  condition  is  sat- 
isfied, both  the  high-frequency  daughter  field  and  the  ion- 
acoustic  wave  can  be  amplified.  The  frequency  matching  con- 
dition for  all  parametric  interactions  which  depend  on  the 
lowest  power  of  the  pump  is : 

These  lowest  order  interactions  couple  the  pumo  wave  to 
two  other  waves  and  are  called  three-wave  interactions. 
Table  I lists  the  possible  three-wave  couplings  in  an  homo- 
geneous, unmagnetized  plasma. 

In  addition  to  these  three-wave  processes,  the  pump 
can  couple  modes  of  the  plasma  throuqh  interactions  in  which 
two  pump  waves  are  involved.  These  processes  are  frequently 
called  second-order  in  the  pump,  even  though  the  phvsics  of 
the  interaction  is  similar  to  the  first-order  processes 
which  lead  to  the  three-wave  interactions.  These  "four- 
wave"  interactions  couple  together  two  pump  waves  and  two 
high  frequency  modes  of  the  plasma  through  a low  frequency 
response  of  the  plasma.  This  low  frequency  response  of  the 
nlasma  has  the  effect  of  making  the  expected  freauency 
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TABLE  I 

THREE-WAVE  PARAMETRIC  INSTABILITIES 


Pump 

Daughter 

Name 

Diagram 

Wave 

Lf  L 


L , A 


L.  T 


T,  A 


T,  T 


2<W 


to0«*2cj, 


Electron- 
ion  decay8 


Uo^UfuL 


C-Ja  Cv* 


Stimulated 
Raman 
scattering 


Stimulated 

Brillouin 

scattering 


11 


Degenerate 
electro- 
magnetic ^ 


matchinq  condition: 


2 (*J piu*p  ^ I €r  -om-  { + | Ld totter  - tn,  I ■ 

hold  only  approximately.  In  contrast  to  the  three-wave 
case,  the  greatest  amplification  docs  not  occur  when  this 
condition  is  satisfied.  These  four-wave  interactions 
are  tabulated  in  Table  II.  In  both  Tables  I and  II,  the 
diaqrams  that  are  used  to  schematically  represent  these 
interactions  are  taken  from  an  analogy  with  quantum  mecha- 
nics. Thus,  the  frequency  matchinq  condition  becomes  a 
condition  quaranteeinq  energy  conservation  in  the  wave  in- 
teraction. 

The  general  formalism  for  parametric  interactions  in 
a plasma  is  well-known  for  the  case  that  the  pump  field  is 
homogeneous,  i.e.,  has  a wavevector,  1^,  equal  to  zero.  In 
the  following  two  chapters,  we  present  a generalization  of 
a kinetic  formalism  of  parametric  instabilities  to  the  case 
o f pumps  with  more  complicated  dependence  on  the  wavenumber 
of  the  pump.  As  we  will  show,  the  most  striking  of  the 
results  of  this  generalization  come  from  the  wavevector 
matching  condition  that  must  be  satisfied: 

fa  A . L 

- o - daughter  - one  ~ tldouqhHr- two  ' 

In  our  quantum  analogy,  this  is  a requirement  for  conser- 
vation of  momentum.  In  general,  this  results  in  important 
changes  in  the  character  of  the  parametric  instabilities, 
which  now  have  a sensitive  dependence  on  the  vector  char- 
acter of  the  pump  wavenumber.  We  have  restricted  our  study 


TABLE  II 

FOUR-WAVE  PARAMETRIC  INSTABILITIES 


Pump 

Daughter 

Name 

Diagram 

Waves 

Lf  L 


T,  T 


T,  L 


L,  L 


L,  T 


Oscillating 


12 


two-stream  ■* 
Filamentation3 3 


Modulational 


U.  ~ U, 


' w,- <•-  =rt/_ 




, - 1-  - ‘k'r 

^ ~ fT 


• ~ <-A 

hJ.Z  Ur 

1 U,  *jj  * 


JO,  =:  Ur 


I C~>I-U)  * Ujr 
t^V+Ul  * l«/i 


Cj,7cUr  5 


L 


T,  T 


-7’ 
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of  arbitrary  k^.,  effects  to  two  general  types:  the  effects 

of  a standing-wave  pump,  and  the  three-dimensional  aspects 
(in  a cylindrical  symmetry  about  kQ) . As  will  be  shown, 
both  of  these  situations  are  frequently  applicable  in  astro- 
physical  and  laboratory  plasmas  in  which  parametric  in- 
stabilities play  a role. 

In  Chapter  II,  the  general  formalism  for  a kinetic 
theory  analysis  of  parametric  mode  coupling  is  derived  for 
a pumo  of  more  creneral  form:  the  pump  can  be  a standina 

wave  (or  travelling-wave)  with  arfcitrarv  wavcvector,  k . 

~o 

This  more  creneral  pump  introduces  two  aeneral  kinds  of 
effects : 

1) .  Effects  due  to  qualitative  changes  in  the  mode 
couplings  possible  with  a standinq  wave  pump,  and 

2)  . Three  dimensional  effects  due  to  the  finite 
size  of  k^j  and  its  vector  character  (independent  of  the 
standinq-wave  character  of  the  pump) . 

The  problem  of  a standinq-wave  pump  driving  paramet- 
ric instabilities  arises  in  several  important  plasma  config- 
urations; whenever  a transversely  polarized  pump  wave  is 
reflected  from  a critical  surface  in  a plasma,  for  example, 
a standing -wave  pump  field  is  set  up  by  the  superposition 
of  the  incident  and  reflected  waves.  A longitudinally 
polarized  standing-wave  pump  field  can  arise  from  the 
decay  products  of  another  parametric  instability  if  these 
daughter  products  have  almost  equal  and  opposite  wavenumbers, 
and  equal  frequencies.  This  is  sometimes  the  case  for  a 
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four-wave  instability  known  as  the  oscillatinn  two-stream 
(OTS)  instability  (see  Table  II).  At  the  end  of  Chapter  II 
we  apnlv  our  creneral  formalism  to  the  specific  case  of  the 
2wpc  and  electron-ion  decay  (EID)  instabilities  driven  by 
a standinq-vave  pump.  The  most  important  effect  here  is 
a lowcrinq  of  the  effective  power  of  the  numo  when  the 
pumn  is  a s tandina-wave . This  reduction  is  a result  of 
different  modes  of  the  plasma  beina  excited  by  each  of  the 
tv.'o  components  of  the  standing-wave.  We  find,  in  addition, 
the  these  modes  are  coupled  more  stronqly  (raisi.na  the 
effective  pumn  power)  as  the  separation  in  freouency  of 
these  waves  decreases  (with  smaller  kQ) , relative  to  the 
arowth  rate  cf  the  instability.  As  kD  qoes  to  zero,  and  a 
standina-vave  becomes  indistinguishable  from  a travellina- 
w?ve,  we  are  able  to  recover  the  travel lino -wave  results. 
When  we  apDiv  our  results  for  the  EID  instability  to 
ionospheric  modification  data,  we  find  the  decrease  in  the 
effective  oumn  power  with  a standinq  wave  nump 
lowers  the  theoretical  predictions  of  the  Lanqmuir  turbu- 
lence produced  in  these  experiments. 

After  this  application  of  our  qeneral  formalism, 
the  third  chanter  is  devoted  to  the  more  complicated,  and 
more  important,  problem  of  the  parametric  instabilities 
driven  by  a longitudinally  nolarized  pump.  The  dispersion 
relation  we  have  derived  makes  possible  a unified  treat- 
ment of  all  Langmuir  pumped  instabilities  which  have  longi- 
tudinally polarized  daughter  waves.  We  find  that  the 


parametric  processes  in  this  case  arc  of  three  tyres,  ill 
of  which  demand  a treatment  in  three  dimensions  with 
arbitrary  wavenumber  pump: 

1) .  /-  non-osci  1 latory  process,  which  wo  have  called 

the  oscillating  two-stream  instability  after  Nishikawa . - " 

Tor  non-zero  k , the  OTS  instability  is  one  dimensional  only 
under  quite  restrictive  conditions,  and,  in  many  cases,  nas 
tts  fastest  orowina  daughter  waves  at  a large  anale  to  Y.  . 

2)  . The  well-known  electron-ion  decay  (LID)  in- 
stability, Tor  k on  the  order  of  the  mass  ratio,  wc  cind 

o 

that  the  LID  instability  can  be  a forward  scattering  nro- 


3).  A new  instability,  which  is  only  seen  in  three 
dimensions,  called  the  stimulated  modulational  (S.M)  in- 
stability. This  instability  produces  dauqhter  waves  at  an 
anqle  to  kQ,  which  can  have  wavenumbers  slightly  larger 
than  kD. 

The  third  chpater  concludes  with  an  apd  ication  of 
these  results  to  three  physical  situations  in  which  the 
longitudinal  pump  is  modelled  on  the  quasilinear  decay  of 
an  electron  beam.  In  all  three  of  the  situations  examined. 
Type  III  solar  radio  bursts,  auroral  arcs,  and  non-relati- 
vistic  laboratory  electron  beams,  we  have  concluded  that 
parametric  effects  do  not  substantially  affect  the  bean 
relaxation,  which  still  proceeds  quasilincarly . other 
aspects  of  these  situations,  however,  do  require  the  para- 


metric processes  of  a Langmuir  pump  for  their  understandina : 


12 


ic  elcctromaanetic  radiation  associated  with  solar  bursts, 
the  anomalous  resistivity  encountered  by  the  auroral  streams, 
and  the  low  frequency  spectra  observed  in  laboratory 
electron  experiments. 

Notation 

We  have  attempted  to  follow  conventional  notation 
whenever  possible.  All  wavenumbers  and  wavevectors  are 
represented  by  lower  case,  Latin  k's;  all  frequencies  are 
denoted  by  lower  case,  Greek  omega.  In  every  case,  a sub- 
script zero  refers  to  the  pump  (for  example,  kQ  is  the  pump 
wavenumber) . We  use  a system  of  units  in  which  time  is 
given  in  terms  of  the  inverse  plasma  frequency  and  distance 
in  units  of  the  inverse  Debye  wavenumber,  kDe,  defined 
above  (note  that  this  gives  velocity  units  of  the  electron 
thermal  velocity).  Finally,  ne  always  is  the  ambient 
electron  number  density  per  cm^  and  Qs  the  temperature,  in 
ergs,  of  the  s species.  Vectors  are  denoted  by  underscoring, 
and  unit  vectors  by  a circumflex;  thus,  the  electric  field 

of  the  pump  is  and  has  direction  c0  or  e (this  is  not 

kQ 

to  be  confused  with  an  unsubscripted  e,  standing  for  the 


electronic  charge) . 
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CHAPTER  II 


THE  GENERAL  FORMALISM 


The  first  studies  of  parametric  instabilities  were 

• 1 2 

done  using  three  simplifying  assumptions:  ' 

1) .  The  pump  wavenumber  is  much  less  than  any  other 
wavenumber  in  the  plasma; 

2)  . A one-dimensional  treatment  is  sufficient  to 
determine  the  fastest  growing  modes;  and, 

3)  . The  pump  is  a monochromatic,  travelling  wave 

field. 


Assumption  one  sometimes  implies  that  the  pump  is 
homogeneous,  the  so-called  "dipole  approximation,"  and 
frequently  is  applicable  to  parametric  instabilities.  From 
the  frequency  matching  condition  for  the  EID  instability, 
(see  Table  I,  page  6),  to  take  one  example,  we  can  derive 
the  following  relation  between  kQ  and  k,  the  daughter 
wavenumber,  for  a transversely  polarized  pumo: 

CkoZ  ^3Vt'k* 


so  that  kQ « k.  If  the  pump  is  a travelling  wave,  this 
allows  us  to  ignore  kQ.  (as  we  shall  see,  assumptions  one 
and  three  are  related,  and  the  condition  kQ « k is  not 
sufficient  to  guarantee  a dipole  Dump  if  the  pump  is  a 
standing-wave) . In  the  derivation  of  a kinetic  theory  of 


parametric  instabilities  in  this  chapter,  we  have  relaxed 
this  assumDtion.  The  necessity  of  studying  more  general 
pump  field  is  clear  if  we  repeat  the  above  calculation 
for  a longitudinally  polarized  Dump,  a Lanamuir  wave.  In 
this  case  we  have: 

k = 2 k-k,t  - ZCt/3  , 

where  cs  is  the  sound  speed  in  the  plasma.  Clearly,  k and 
kQ  can  be  comparable  in  maonitude.  We  have  found  the  most 
important  effects  of  a finite  kQ  pump  do  occur  for  lonai- 
tudinally  polarized  pumps.  In  Chapter  III  we  v*ill  derive 
detailed  results  for  the  kQ  dependent  effects  in  para- 
metric effects  pumped  by  a Langmuir  wave  field. 

Furthermore,  when  kQ « k,  we  would  expect  that  the 
parametric  processes  would  not  depend  on  the  direction  of 
kQ  (except  as  this  is  determined  by  the  polarization  of 
the  pump),  and  hence,  that  a one  dimensional  treatment 
would  be  sufficient.  This  is  the  content  of  assumption 
two.  However,  to  be  consistent,  if  we  relax  assumption  one, 
we  must  be  ready  to  examine  possible  three  dimensional 
effects.  Indeed,  our  formalism,  especially  in  the  reaimes 
when  kQ  ^ k,  predicts  important  three  dimensional  effects 
for  all  parametric  instabilities.  These  will  be  examined 
in  detail  in  Chapter  III,  for  the  case  of  a Langmuir  pump. 
Tot  the  EID  instability,  to  compare  with  our  example 
above,  driven  by  a Langmuir  pump,  not  only  are  k and  kQ 

comparable,  but,  in  addition,  the  anqle  between  k and  k 

— 


i 


1C 

enters  critically  into  the  condition  for  maximum  growth. 
Our  general  formalism  is  derived  in  a three  dimensional 


| 


geometry  in  which  we  have  assumed  a cylindrical  symmetry 
about  Kq. 

Finallv,  our  motivation  for  relaxing  the  assumption 
concerning  a travelling  wave  pump  field  (assumption  three-- 
although  we  have  kept  the  assumption  of  a monochromatic 
Dump)  is  the  evidence  that  in  many  experimental  situations 
the  pumo  is  not  a simple  travelling-wave.  It  is  frequently 
the  case,  for  examole,  that  the  theory  of  parametric  in- 
stabilities is  apolied  to  experiments  in  an  inhomoqeneous 
plasma.  In  this  situation,  one  approximation  to  a treatment 
of  the  problem  which  take  the  density  variation  into 
account,  is  based  on  the  observation  that  an  electromag- 
netic pump  will  be  reflected  in  the  region  of  anv 
critical  density  surfaces  in  the  plasma  (where  the  local 
plasma  frequency  equals  the  frequency  of  the  pump  wave) , 
and,  the  superposition  of  the  reflected  and  incident  waves 
will  create  a standing-wave  electromagnetic  field  in  thr 
plasma.  If  the  absorption  of  energy  from  the  pump  is  small 
enough,  this  standing-wave  will  double  the  local  amplitude 
of  the  pump  field,  but  with  a new  wavenumber  and  frequency 
dependence.  This  is  thought  to  be  the  case  for  the  plasma- 
radiation  systems  occurring  in  ionospheric  modification  ex- 
periments and  in  laser-irradiated  pellets. 3 The  other 
general  class  of  problems  in  which  a standing-wave  pump 


.j 
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arises  is  that  of  the  saturation  of  parametric  instabil- 
ities. One  example  of  this,  which  will  be  examined  in  more 
detail  in  the  next  chapter,  is  the  non-oscillatory 
interaction  of  a small  wavenumber  pump  whose  daughter 
waves  are  a pair  of  Langmuir-like  waves  with  almost  equal 
frequencies  and  equal,  but  opposite,  wavenumbers. 

The  superposition  of  these  two  dauqhter  waves  sets  ud  a 
standinq- wave  field,  and  the  saturation  of  the  process 
depends  then  on  the  dynamics  of  the  standinq-wave  pump 
generated  by  the  oriqinal  dauqhter  waves. 

In  this  chapter,  we  beqin  with  a kinetic  theory 
derivation  cf  parametric  excitation  in  an  homogeneous, 
infinite  maqnetic  field-free  plasma.  Our  derivation 
relaxes  assumptions  one,  two  and  three.  We  have  divided 
our  applications  of  this  formalism  into  two  parts: 

1).  At  the  end  of  this  chapter,  we  study  the  effects 
denendent  on  the  standina-wave  pump.  The  most  strikino 
effect  we  find  is  a lowering  of  the  effective  power  of  the 
pump.  This  is  an  effect  we  explain  by  a detailed  consider- 
ation of  the  new  mode  couplings  that  a standinq-wave  allows. 
These  new  couplinqs  arise  from  the  different  wavenumber 
matching  conditions  possible  with  each  component  of  the 
standinq-wave.  The  key  concept  concerning  these  phenomena 
is  that  of  the  "closeness"  of  modes;  that  is,  the  separation 
in  frequency  of  two  modes  compared  to  the  growth  rate  of 
the  instability.  tisinq  this  idea,  we  have  derived  rela- 
tions between  kQ,  the  number  of  modes  coupled  in  the  plasma. 


and  the  effective  power  of  the  nump  These  results  are 
applied  to  ionospheric  modification  experiments,  where  our 
predictions  of  a lowered  effective  pump  rower 
changes  the  theoretical  predictions  made  of  the  observed 
Lanqmuir  turbulence. 

2) . All  of  Chanter  III  is  devoted  to  a study 

the  three  dimensional  and  finite  k effects  of  parametric 

o 

instabilities  driven  by  a Lanqmuir  wave  pump. 

THE  KOPsMAL  MODE  EQUATION4 

Maxwell's  equations  in  a medium  are  efficiently 
solved  beginnina  with  the  equation  relatino  the  ensemble- 
averaoed  electric  field,  F,  the  external  currents,  and  that 
part  of  the  internal  current  which  has  a hiqher  than  linear 
dependence  on  the  field,  ; 

(1)  *}(ktu>)E(k,v) ♦ j"L(k, (*))], 

where  3 (k,  ) is  the  generalized  susceptibility,  related 

to  the  linear  dielectric  function,  C(k,  w)  through  the 
relation : 

wm  / 1 . — fern  c1  fe*  . r / 1 v 

(2)  "l/m  (>!,“) £ UCtmiM. 

The  dielectric  function  itself  is  derivable  from  the 
equations  describing  the  plasma.  Our  strategy  in  solvinq 
this  equation  is  as  follows: 


1).  Ev  limiting  ourselves  to  linear  stability  ana- 
lysis, to  look  at  the  solutions  to  the  equation  when  there 
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are  no  external  currents,  jext  = Qt  and  concern  ourselves 
only  with  the  normal  modes  of  the  system,  whose  dispersion 
relation  is  given  by  roots  of  : 

(3)  12  (M  £(*,“)  = -4ft.'  j 

2) .  To  find  an  expansion  of  the  non-linear  current 
in  powers  of  the  field;  and, 

3)  . To  introduce  an  orderina  among  the  contributions 
to  the  field,  distinguishing  between  the  pumn  field,  an 
externally  applied  field  which  is  assumed  to  be  of  fixed 
amplitude,  frequency,  and  wavenumber,  and  the  much  smaller 
fields  that  the  pumo  drives  up  (this  is  the  so-called  para- 
metric approximation) . 

The  solution  to  (3) , using  only  the  internal  current 
which  is  first  order  in  the  pump,  i.e.,  second-order  in 
the  total  field,  gives  the  three-wave  parametric  interac- 
tions; the  current  which  is  third-order  in  the  total  field 
(second  order  in  the  pumo)  gives  rise  to  the  four-wave 
processes.  To  solve  equation  (3),  we  must  evaluate  the 
internal  current;  in  general,  we  can  write  the  current  in 
terms  of  the  non-linear  susceptibilities  and  the  fields  in 
the  plasma.  The  linear  susceptibility  is  given  by  the 
linear  dielectric  function; 

(4)  j/°  (!Siw)  * ~ (k,v)  E 

where : 

= (&>•*>)  • 
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The  non-linear  currents  can  likewise  be  written  in  a 
series  expansion  which  comes  out  of  a Fourier  transform  of 
an  expansion  of  the  current  in  terms  of  the  electric  field: 

j*(M  * ju>(M+  + ••• 

where : 

(5)  O*!**1)  “ * ^ (£,w,  ltu/x)  E»(kz,Ui) , 


and : 


(6) 


jl*(ktu)  - ~lu>  j dKin  fanp  ki.Vi,  *»#•»*) 

* En(^i,0Jt)  fUli) 


with  the  definitions: 


7)  JKn  ■ ■--fy?*'  JU‘  S'Os-k-k)  %(u-U.-ut), 


(8) 


i k.«  » sn-k-k-k)  »(ww. 


It  is  important  to  note  that  the  model  of  the  plasma-field 
interactions  used  to  derive  the  susceptibilities  is  sepa- 
rate from  the  model  used  to  evaluate  the  dielectric  function 
in  equation  (2).  We  have,  in  the  Appendix,  evaluated  JLfa* 
and  Xfa*  for  several  important  cases  using  a fluid  model.4 
From  equations  (A. 7)  and  (A. 8)  in  the  ADpendix,  we  have 
explicit  forms  for  the  non-linear  susceptibilities,  and 
we  use  them  with  the  definition: 


(9)  7 Cktui  4'jW,  j |?t ,UJi]  * e*C*,en  Xfm*  > 

where  cj  = e^,  em  = ekl,  etc.  are  the  unit  polarization 
vectors  of  the  E(k) , E(k^) , etc.  The  next  step  is  to  make 


J 


and  that  EQ  still  be  small  enoucrh  to  ensure  the  convergence 
of  the  current  expansion  (a  sufficient  condition  for  this 
convergence  is  that: 

f |E.|*  d'k  « »e£  . 

whose  details  are  discussed  in  reference  5) . Physically, 
our  zero  order  (in  E')  state  is  that  of  a large  amplitude 
known  field  interacting  with  a plasma.  Our  object  is  to 
calculate  the  effects  linear  in  E'.  In  general,  tQ  has 
the  effect  of  coupling  together  previously  independent 
modes  of  the  plasma  (resonances  in  E')  and,  in  some  cases 
amplifying  these  now  coupled  resonances.  This  is  a para- 
metric process  since  we  take  |ec|2  as  fixed,  and  the  time 
and  space  variation  of  EQ  provide  the  coupling  mechanism. 

We  now  relax  the  assumptions  listed  in  the  beginning 
of  this  chapter,  on  the  form  of  the  pump,  and  write: 


(11) 
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where  f is  an  arbitrary  phase  factor.  This  pump  is  a standir.g- 
wave  of  arbitrary  kQ  when  Ot  ==  S and  a travelling-wave  when 
K = 1 or  0,  The  local  amplitude  of  this  pump  is  always  J E0J 
for  any  value  of  0(.  We  write  the  Fourier  transform  in 


space  and  time  of  the  pump  as : 


Substitution  of  this  expression  for  E0  in  (10)  and  then 
into  (5),  (€),  and  their  higher  order  analogues,  immediately 
gives  the  various  currents  in  (4) . The  second  order 
current  is,  for  example: 

+ ♦ p. ^ (fe , W) * 3 jl., 5 k,,  to • ^ 

In  an  isotropic  plasma  (BQ  = 0)  , we  can  choose  m(k,w) 
diagonal,  and  equation  (3)  becomes: 

d3)  u)E/ = - ~j~  [ji  (%>&)  + 

This  equation  generates  a system  of  equations  connecting 
fields  at  various  values  of  the  wavenumber  and  frequency, 
shifted  in  their  arguments  by  letting  k go  into  kikQ  and 
wgo  into  0JtuQ.  The  result  is  an  infinite,  homogenoeus 
set  of  coupled  equations  for  E’ (k±nk0,  <JiwQ)  . 

Physically,  the  pump  has  coupled  previously  indepen- 
dent modes  of  the  plasma;  there  is,  for  example,  the 
term  in  these  equations  proportional  to  which  couples 

the  mode  at  (k,w)  with  the  four  modes  at  (kikQ,  0>ia/o) 
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and  w?u0)  . Our  object  is  to  solve  this  system 

of  equations.  Given  that  we  know  the  susceptibilities, 
substitution  of  (12)  and  the  like  expressions  for  etc. 

into  (13)  qives  a homoqeneous  set  of  equations  for  the 
b'  fields  at  different  wavenumber  and  frequency,  and  the 
determinant  of  the  coefficient  matrix  of  this  set  of 
equations  qives  us  a dispersion  relation  which  describes 
the  n ew  normal  modes  of  the  plasma.  We  now  turn  to  an 
evaluation  of  the  dispersion  relation. 

THE  DISPERSION  RELATION 

The  fundamental  difficulty  in  solvino  the  matrix 
equation  that  follows  from  (13)  is  that  it  is  infinite. 

The  infinite  number  of  variables  that  we  have  are  the 
fields  E':  E'  (k,  U)  , E’fkikQ,  W±fc>Q)  , E ' (k±k0,  w»«0)  , 

E'  (h±2k0,  oJ) , etc.  We  must  find  a way  to  truncate  this 
doubly  infinite  qroup  of  variables  in  wavenumber  and 
frequency,  if  we  are  to  get  a finite  set  of  equations.  We 
will  show  that  truncation  in  frequency  and  wavenumber  is 
provided  by  a restriction  that  allows  us  to  take  into  account 
only  modes  that  fall  near  a resonance  of  the  unperturbed 
system.  However,  our  ability  to  truncate  the  series  of 
variables  in  wavenumber  will  be  a sensitive  function  of 
)£0,  which  is  not  the  case  for  the  truncation  in  frequency. 

In  this  section  we  will  reduce  our  infinite  equation  to 
an  m x m determinant,  where  m is  a function  of  k0.  We 


then  will  apply  this  m x m determinental  dispersion 
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relation  to  several  parametric  instabilities  driven  by  a 
standina-wave  pumo. 

We  now  turn  to  the  problem  of  truncatinq  the  infinite 
matrix  equation  that  follows  from  (13).  We  first  make  the 
convention  that  0 < w <w0.  This  guarantees  that  all  the 
arguments  in  frequency,  utnUg  are  different  for  different 
values  of  n.  One  consequence  of  the  condition  that 

f I Eol  « n Ge  , 

is  that  the  predominate  response  of  the  plasma,  E',  will 
be  very  near  the  normal  modes  of  the  unperturbed  plasma. 

It  is  in  this  sense  that  the  pump  is  said  to  couple  the 
normal  modes  of  a plasma;  that  is,  if  W,  is  a high  fre- 

a 

quency  normal  mode  of  the  plasma,  so  that: 

€(K»  i ^ 0, 

and,  if  w(k)  is  a solution  to  our  new  dispersion  relation, 
then : 

)w  - £da|  <<  6d0  • 

This  is  easily  verified  a posteriori.  Now  the  normal  modes 
of  the  unperturbed  plasma  are  well  known,  and  we  can  con- 
clude that  we  need  only  keep  fields,  E',  which  have  a 
frequency  near  a normal  mode  of  the  plasma.  Since  we  will 
choose  the  pump  frequency,  wQ,  to  be  a high  frequency  (on 
the  order  of  the  plasma  frequency) , it  is  easy  to  verify 
that,  for  the  responses  of  a homogeneous  plasma,  Table  I 
(paqe  6)  gives  the  possible  excited  fields  for  the  couplinas 
allowed  by  j^).  Let  us  call,  for  such  a three-  wave 


i 
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process,  wa(k^)  and  ^(kj-,)  - • the  normal  modes 

of  the  unpumped  plasma;  the  dauqhter  waves  of  the  insta- 
bility must  be  near  these  two  frequencies.  Thus,  for 
example,  to  excite  the  2wne  instability  (see  Table  I),  we 
must  satisfy  the  frequency  matching  condition,  at  least 
approximately : 

U)  S:  6 Jl  (.ka)  S U>a  ClSa)  > 
u>  - cj„  ^ - u/L  (A*  -(?«)  = t*'1*  (k*-k o) 

(in  the  case  of  the  2WpG  instability,  u>a  and  the  natural 

resonances  of  the  plasma  that  the  pump  excites,  are  both 
Langmuir  waves) . We  can  choose  Jia  so  that  these  equations 
are  satisfied.  But,  the  other  modes,  at: 

U)  + UJ0  ~ 3d/|u_, 
co  - 2.W,  ^ ~2 

will  never  be  near  a resonance  of  the  unpumped  plasma.  The 
more  exact  meaning  of  "near  a resonance"  is  straioht- 
forward.  Physically,  we  want  w(k),  the  new  normal  mode,  to 
be  within  a growth  rate  of  the  "old"  normal  mode  of  the 
system.  Expressing  this  in  terms  of  the  dielectric  fun- 
ctions, which  are  more  properly  a measure  of  the  plasma 
response : 

|R*eu,u>)-R*eck,u,,)|  < | e (&,<*>)!  , 


Fiqure  one  a)  shows  schematically  what  is  required.  These 
inequalities  allow  us  to  truncate  our  equation  in  frequency. 


26 


Figure  one.  Qualitative  diagram  of  the  plasma 
response  as  a function  of  wavenumber  and  frequency.  In  a) 
we  have  the  criterion  for  determination  of  which  frequencies 
are  "close  enough"  to  a natural  resonance  of  the  plasma 
.(shown  by  the  dotted  lines  at  U (ka)  and  -^(kfc,)).  The 
mode  at  *J-«0  is  close  enough;  the  mode  at  0>-2fc>o  is  not. 

In  b)  the  criterion  for  determination  of  which  wavenumbers 
will  be  excited.  Modes  at  (ka±kQ,  u-u>0),  (ka,w),  and 
(ka-2k  , u)  are  "close  enough"  to  a natural  resonance,  but 
the  modes  at  (ka-3k0,  w-uQ)  and  (ka+2k0,  u>)  are  not. 
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Notice  that  if  w is  a low  frequency,  that  it  is 
possible  to  have  three  daughter  frequencies  near  resonance: 

£d ' l*/«  ~ - Df*.  } CjO+  « LOfkA.  , CO  « U)A  . 

This  enables  E'(k,  w)  , E'(k+ko,  £o  + wQ)  and  E'(k-k0,  <o-wQ) 
all  to  be  excited.  We  shall  see  that  is  the  case  for  the 
four-wave  instabilities  (see  Table  II,  page  8).  in  the 
remainder  of  this  chapter,  however,  we  assume: 

|Rx.  - R*  6(cva)|<  kr\t(u/-CJ,)  <•  I & £(u>*Ulo)  ~ |, 

and  we  will  defer  discsssion  of  four-wave  effects  to 
Chapter  III. 

We  now  specialize  our  results  to  this  case,  and  make 
the  related  assumption  that  only  needs  to  be  consider- 

ed. This  is  admissible  as  long  as  we  can  neglect  small 

4 

corrections  to  Re  u,  and  are  interested  only  in  three-wave 
effects.  In  chapter  III  we  will  include  effects  due  to 
both  the  up-shifted  wave  at  E' (k+kQ,  u+uq)  and  due  to 
couplings  introduced  by  j(3).  This  leaves  us  with  the 
following  equations,  obtained  from  a substitution  of  (12) 
into  (13) : 


k 


kitAf)  E„  (k-ko 

' 1 Mr 

k,u)  Bn'  (h.lk*,b-<*>i)Eor*  i 

♦ k'kj>,tJ-io0)  £„' (k-Zk.,u) 

(16)  *»> «*>  ^ E;  (k*2k.,u)EtM 

♦(I-*)  En'(k,u)  i 

(17)  i(X-lk.tu)  Ej  (fe-2k,,fc;)  - <x (Jt- 3 

+(,_^  k*Mi  K-2 ^u)  (k-fe.w-w.)  E0t ; 
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(16)  (k*2ko  ,(S)  — (X  fc^'^-fU/M*,'  k4?kc,i«)l>  (k*ke,^  (A^jE^n. 

4ft 

+ 0 -^/'JLmm.  (k<  IkoUH*;  k*2  k,  w)E,'(f^2fc.,u)  ten 

When  o<  = 0 or  1,  and  we  have  a travelling-wave  nunp, 

only  equations  (14)  and  (15)  are  necessary.  That  is,  with 
a travelling-wave,  a truncation  in  frequency  also  truncates 
the  series  of  equations  in  wavenumber.  This  is  a result 
of  the  fact  that  a travelling-wave  can  only  couple  modes  in 
which  the  wavenumber  and  frequency  are  shifted  by  the  same 
multiple  of  (kQ,  WQ)  . However,  a standing-wave  can  couple 
modes  a more  complicated  fashion.  Equation  (19)  shows,  in 
matrix  form,  the  infinite  set  of  equations  qenerated  like 
(14) -(18)  by  a standinq-wave  pump.  These  equations  have 
been  truncated  in  frequency,  but  not  yet  in  wavenumber.  An 
analogous  truncation  is  possible  for  the  shifts  in  wavenum- 
ber, since,  for  large  enough  kQ,  efkfko,  lj-u>0)  and 
6(k+3kc,  io~ Uq ) cannot  both  be  on  resonance  simultaneously. 
That  is,  if  we  assume  that  £(k-k0,  fc-ofc,)  is  very  close  to  a 
resonance  of  the  unpumped  plasma,  then  the  mode  at  (ktnkQ) 
with  frequency  (t>-u0)  will  be  "near"  this  same  resonance 
only  if: 

|R*  ecfc.w)  -R£e(k*r)ak,u;)|  < I 3**  e(fe,w)| » 

or,  in  terms  of  the  other  mode: 

< ItV* c(k-kc, u>-u>.)\  . 

Our  wavenumber  and  frequency  matching  conditions  and 
eouation  (19)  allow  these  equations  to  be  written: 

(20a)  KM  IStf/Wa]  | < |JU  6 (k,  | 


A 


(20b)  ) Rt  e(ka 'h/u\)-Ra.c(&*:nfa>i‘'J*)l  < I Jva  efk-^W'Wc)! . 


notice  that  the  same  integer,  n,  enters  both  (20a)  and  (20b). 
This  reflects  a physical  property  of  (19),  namely,  that  two 
conditions  must  be  satisfied  for  a given  mode  to  be  ex- 
cited (thus,  not  excluded  by  our  truncation  scheme): 

1) .  The  ride  must  be  "near  enouqh"  to  a resonance 
of  the  unpumped  system  (this  is  expressed  in  (20a));  and, 

2) .  This  mode  must  be  coupled,  in  (19),  to  another 
mode  that  satisfies  the  first  condition. 

This  second  requirement  is  a result  of  the  fact  that  the 
second-order  susceptibilities  can  only  couple  modes  shifted 
by  (±kQ,  and  both  of  these  modes  must  be  close  to 

resonance.  For  example,  in  the  EID  instability  (see  Table 
I,  page  6),  when  u>a  is  a Langmuir  mode  and  ^ is  an  acoustic 
mode,  the  na  and  nb  may  be  different.  We  must  pick  the 
smaller  of  na  and  nb,  i.  e.,  the  one  that  satisfies  (20a) 
and  (20b)  simultaneously  (na  and  nb  are  defined  in  the 
two  inequalities  immediately  preceeding  (20a)  ar.d  (20b)  ). 

In  the  next  section,  we  will  show  results  from  a numerical 
test  of  this  second  requirement,  by  evaluatinq  (19)  in 
the  case  when  we  use  (20a)  and  (20b)  to  eliminate  modes,  and 
in  the  case  when  we  evaluate  (19)  without  taking  the 
smaller  of  na  and  n^. 

equations  (20a)  and  (20b) , then  allow  us  to  truncate 
(19)  ir.  wavenumber.  Figure  one  b)  shows  the  situation 
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where  modes  at  (katkQ,  u-o0)  , ():a,W),  and  (ka-2k0,«M 

will  be  coupled  together.  The  other  modes  shown,  at 
(ka+2kQ,  c*>)  , (ka-3k0,  V-uQ)  are  not  "near  enouqh"  to  a 
natural  resonance  of  the  system  to  be  excited.  In  most 
problems  of  interest,  all  frequencies  of  the  daughter  waves 
will  be  functions  of  the  wavenumber  which  are  slowly  varyino, 
and  so  we  can  write  (20a)  and  (20b)  as: 


(21a) 


(n-i)  foe 


36  ((;,') UJ) 

" 1*  k > 


£ (&*,*>)  7 


Jkiw 6 Lk-'k*' ~lct J I • 


(21b)  ' 3k' 

where,  again,  n is  the  largest  positive  integer  for  which 
both  inequalities  hold.  Examining  (21)  or  (20)  we  can 
determine  the  total  number  of  modes  that  a pump  with  a 
given  kQ  can  couple  together: 

1) .  If  n ■ 1,  and  |k'a  (ka)  +^b(^a-ko^  ~ ^o  ^ 

we  cannot  satisfy  the  basic  frequency  matching  condition. 

In  the  three-wave  approximation,  there  are  no  parametric 
interactions  possible  for  this  pump  (note  that  this  case 
could  also  arise  if  there  were  no  value  of  n,  besides 
zero  perhaps,  which  could  satisfy  (20a)  and  (20b)  together). 

2) .  If  n = 1,  and  I t/a(ka)  +Wj3(ka-k0)  - |<  Jlm^j, 

we  have  two  possibilities.  If  equation  (20b)  holds  only  for 
the  upper  sign,  then  the  following  modes  are  coupled: 

4 blk),  • 


m— — — —ft—i  I— ma— I— 
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This  is  the  sane  couplinas  that  exist  in  the  usual  travel- 
ling - wave  pumped  instability.  If,  however,  (20b)  is  satis- 
fied for  both  signs,  then  the  following  modes  are  coupled: 

fod(k)  Vb{Mi60)- 

This  is  the  simplest  standing-wave  coupling. 

3).  If  n>l  (for  both  signs  of  (20b)),  then  we 
couple  together  m modes,  where  m = 2n  + 1 . 

Thus  the  inequalities  (20)  or  (21)  allow  us  to 
truncate  the  infinite  matrix,  (19)  in  wavenumber,  keeping 
only  an  m x m matrix  out  of  the  center  (the  2x2  matrix  and 
3x3  matrix  described  by  case  two  above  are  marked  by 

dotted  lines  on  equation  (19)).  In  the  2x2  case,  we  get 
immediately  the  determinental  condition: 

(22a) 

- £ k,fco-w0)/  (k'k*,ui-vt-£e;u>,;ktv)  |E*|7  r0. 

When <<  = 1,  this  is  the  result  of  the  usual  travelling-wave 
calculation.  Notice  that  for  a standing-wave  pump,  when 
* = *5,  we  have  a lower  effective  pump  field.  The  conse- 
quences of  this  will  be  examined  in  the  next  section.  For 
the  3x3  determinant: 

This  is  the  dispersion  relation  we  will  study  in  detail  in 
the  next  section,  applying  it  to  the  cases  of  the  EID  and 
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2t)pe  instabilities.  One  feature  of  this  equation,  which 
bears  on  the  meaning  of  our  truncation  scheme,  should  be 
noted;  (22b)  implies  (22a)  whenever  (20b)  is  not  satis- 
fied for  both  signs.  That  is,  if  £( k+is0»  U>-u0)»  £ (k-kQ,  v-uQ) 
which  is  equivalent  to  (20b)  being  violated  for  the  upper 
sign,  then  (22b)  simplifies  to  (22a).  This  is  a result  of 
the  fact  that  our  truncation  scheme  does  not  introduce  any 
new  physics;  it  only  serves  to  make  the  problem  mathemati- 
cally tractable  by  giving  us  a method  of  predicting  which 
modes  will  be  important  (i.e.,  coupled  toqether)  and  which 
we  can  afford  to  ignore. 


APPLICATIONS 

We  now  apply  our  general  formalism  to  specific 
situations  in  which  a transverse  pump  decays  into  longitu- 
dinal waves.  There  are  two  possible  three-wave  processes 
in  which  this  happens: 

1) .  The  2 c»pe  instability.  Here  a pump  wave  near 
twice  the  plasma  frequency  decays  into  two  Langmuir  waves. 
We  will  show  that  the  geometry  and  growth  rates  of  this 
instability  are  changed  considerably  by  a standing-wave 
pump.  The  qualitative  changes  are  due  to  the  superposition 
of  two  instabilities,  each  due  to  one  component  of  the 
standing  wave. 

2) .  The  EID  instability.  In  this  instability,  the 
pump  wave  decays  into  an  acoustic  wave  and  a Langmuir 
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wave.  As  was  noted  before,  this  instability  always  has 

k >>  k0  and  the  effects  of  a standing  wave  are  confined  to 

a decrease  in  the  effective  power  of  the  pump.  Even  so, 

an  analysis  of  the  EID  instability  in  the  interpretation 

of  ionospheric  modification  experiments  shows  this 

effect  to  be  significant.  For  both  instabilities  we  examine 

the  limit  of  small  k to  recover  the  standard  travellinrr- 

o 

wave  results. 

The  2 u pc  Instability^ 

The  2Wpe  instability  can  be  visualized  as  the  decay 
of  the  pump  wave,  here  at  near  twice  the  plasma  frequency, 
into  two  plasma  waves.  The  approximate  freauency  matchina 
condition  for  this  decay  process  is: 

(22)  We(fco)  Idi(Jl)  -t  Ul 

This  matching  condition  and  the  results  of  ineauality  (20) 
determine  the  qualitative  aspects  of  the  instability.  Ir. 
the  initial  work  on  this  instability,  Goldman  derived 
results  for  arbitrary  kQ.6  The  2ope  instability,  in  fact, 
does  not  exist  for  a homogeneous  pump.  He  found  maximum 
qrowth  when  (22)  is  satisfied  and  k'kQ  l/f~2.  The  only 
generalization  our  results  offer  is  that  of  a standir.g- 
wavc  pump.  We  find  two  regimes  in  the  effects  a standing- 
wave  pump  produces.  One  is  the  regime  in  which  the  dauqhter 
waves  pumped  by  each  component  of  the  standing-wave  arc 
independent;  this  results  in  two,  travellinq-wave  in- 
stabilities, pumped,  however,  by  a field  of  only  half  the 
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local  amplitude.  In  this  case  we  recover  the  geometry  of 
the  travelling  wave,  where  k*kc  s=*  1/fT.  In  the  other 
regime,  the  daughter  waves  interact.  This  results  in  a 
growth  rate  intermediate  between  the  other  standing-wave 
regime  and  the  growth  rate  that  would  result  from  a travel- 
ling wave  pump.  The  geometry  is  essentially  unchanged, 
rinally,  we  use  physical  armaments  to  show  that  in  the 
k >>  kc  regime,  the  standing-v/ave  and  travelling-wave  results 
agree.  Let  us  look  at  the  simplest  case  of  a standing- 
wave  drivinq  a 2Wpe  instability,  v/hen  m = 3,  and  our  dis- 
persion relation  is  given  by  (22b).  In  this  case,  we  have 
the  following  three  modes  coupled: 


<23>  foioCWs^t*);  . 


and,  from  (20)  we  get  the  exact  conditions  on  kQ: 

(24)  3k,  • (fcffc)  < < 6ko  (k*2ko)  ■ 

Usinq  the  symmetry  properties  of  the  susceptibilities  given 
in  the  Appendix,  the  equation  describing  the  new  normal 
modes  of  the  plasma  is  (from  (22b)): 

£( fe,w)  e(k -*»,  u-  v.)  €(k*ko,  to-u)  t ex'  (k,v;  k- ko,u,-u)l\ 

(25)  € Ck*hu, 


The  Appendix  gives  the  general  form  for  these  susceptibilities: 


(26) 


3b 


where : 


t«  We,, 

W0  = E'V[4tmc(&c  * 6i)J  . 


Since  all  the  high-frequency  modes  arc-  well-defined,  we 
also  make  resonance  approximations  for  all  the  dielectric 
functions ; 


(27a)  e(k,U>)  - [(0-  UK  (D*  i?f]  , 

(27b)  + • 


This  reduces  equation  (25)  to  an  algebraic  equation,  cubic 
in  the  frequency: 

[w-U}l(k)*i)f][k)-u/e+£>t(k+6»)  + Tx][v-Ho+  Ui  (k-Mc )+>'y] 

(28)  ~OLir.Z[u-CJo+U)L(k+kt)<it'] 

-(l-ot)irS[uJ-w,4  Ui(b-k.)*i)i]  = 0. 

Although  it  is  not  necessary,  we  will  also  assume  that  the 
only  damping  of  these  modes  is  collisional,  and,  so  not 
dependent  on  k.  This  cubic  equation  is  intractable  in  its 
full  generality,  but  it  can  be  solved  approximately  in  two 
regimes  where  the  instability  exists.  Defining  the  fre- 
quency mismatch  factors: 

St  = U)9  - H (k  i kc)  f 

we  have  the  following  conditions  under  which  the  solutions 
to  equation  (28)  will  have  a positive  imaginary  part  and  can 
be  dealt  with  analytically: 

A).  Independently  excited  daughter  waves.  In 


this  case,  one  of  the  daughter  waves  is  much  nearer  to 
being  resonant  than  the  second,  or: 

Jn\  £ > C (k»± k ,00-ujc)  » £ (kc  , to-  u)o) 


I | >>  I (a)  - St  | • 

Equation  (28)  can  be  approximately  written: 

(29)  [w-  uL(k)*  ii[[b)  - S±  + ;'y]  + ^ =0. 

This  is  the  equation  describing  the  usual,  travelling-wave 
pumped  instability,  except  that  the  coupling  between  the 
two  Langmuir  waves,  P+,  is  only  one-quarter  the  coupling 
for  a travelling  wave  pump;  this  is  a result  of  the  local 
amplitude  of  the  effective  pump  being  one-half  EQ.  In 
other  words,  only  one  component  of  the  standing-wave  pump 
drives  the  instability.  The  growth  rate  is  then,  solvinq 
equation  (29) : 

C30)  = -Jf  *J  j [It T j » 


which  has  its  maximum  value  of: 


=f(,;<<)}  r*  -y- 


when  (^(k)  = Notice  that  the  maximum  growth  occurs 

when  equation  (22)  is  satisfied;  this  is  the  condition  that 
determines  k for  maximum  growth,  as  a function  of  k^  and 
so  we  have  deduced  the  frequency  matching  condition  in  a 
natural  way  out  of  the  dispersion  relation.  In  addition, 
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there  is  a regime,  for  values  of  ft  very  near  the  threshold 
value  of  I"*  =4/1  , at  which  the  maximum  growth  rate  is  pro- 
portional to  P*  rather  than  P.  Figures  two  and  three  show 
a numerical  solution  to  equation  (28)  in  the  regime  of 
case  A.  Figure  two  shows  the  real  and  imaginary  parts  of 

A A 

W,  with  k-k0  = 45°,  and  Figure  three  is  a plot  in  polar 
coordinates  of  the  imaginary  part  of  the  frequency  (the 
growth  rate  of  the  instability)  as  a function  of  k.  Notice 
that  for  this  case,  equation  (20),  which  requires  the 
modes  at  kik^j  to  be  "close  enough"  to  resonance,  is  at 
best,  marginally  satisfied;  we  have,  physically,  two,  es- 
sentially uncoupled  (except  near  90°)  instabilities,  each 
pumped  by  a travel ling -wave. 

B) . Coupled  excited  daughter  waves.  When  6+=  S_, 
both  daughter  waves  can  be  near  resonance  for  the  same  k, 
and  equation  (28)  becomes  (for  t*=  1-&)  : 

(3D  [w-u>t(k)  ♦»*■][  w - S+  ij  ] h 2**r  -0 

where : 

$=  (S*+  s.)/z 

r*  (rvr.)/2. 

The  maximum  growth  rate  here  is  larger  by  *{2  than  for  case  A 

(32)  <9rn  wL  « r , u>Uk)  - S . 

* w 

Figure  four  shows  the  real  and  imaginary  parts  of  the  solu- 
tions to  equation  (28)  in  this  regime,  and  Figure  five 
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Figure  two.  Independently  excited  daughter  waves 
driven  by  a standing -wave  pump  at  $>kQ  = l/f2.  The  real 
parts  of  the  three  roots  of  (28),  given  in  the  top  solid 
lines,  and  the  imaginary  part  of  any  root  which  greater 
taan  the  Langmuir  damping,  are  shown  on  the  same  graph. 

The  following  parameters  put  this  example  in  the  regime  of 
case  A of  the  text  for  the  2 pe  instability:  kQ  = .056; 

w0  = 1.6x10-4;  (ve/c)2  = 10-3;  and  < «=  10"3.  The  maxi- 

mum growth  rate  occurs  when  (23)  is  satisfied.  In  this 
regime,  equation  (29)  is  a good  approximation  to  the 
above  solution  in  the  regions  of  growth. 


Figure  three.  Growth  rate  as  a function  of  k for 


independently  excited  daughter  waves  in  the  2 wpe  instability. 
This  is  a polar  contour  plot  of  the  imaginary  part  of  the 
solutions  to  (28).  The  symmetries  come  from  the  mathe- 
matical properties  of  (29),  such  that,  for  example,  in- 
verting k chances  the  upper  for  lower  sign  in  (29).  The 
signs  in  the  center  of  each  region  of  the  instability  cor- 
respond to  the  appropriate  sign  in  (22)  and  also  indicate 
the  wave  coupling  that  produces  the  instability  (i  cor- 
responds to  w(k)  and  u(k±k0)  being  coupled).  The  para- 
meters for  this  graph  are  the  same  as  in  Figure  two.  The 
units  in  the  contours  are: 

0 -lxl0“3  U)pe 

10 

2 +lxl0~3  UpC. 


Figure  four.  Coupled  excited  daughter  waves  driven 
by  a standing-wave  pump  at  lt*K0'/fi,the  2Wpe  instability.  The 
real  roots  of  (29)  are  given  on  top.  and  the  imaginary  part 
of  any  root  when  it  is  positive,  given  in  the  lower  curve. 

The  following  parameters  put  this  example  in  the  regime  of  n 
case  B:  kQ  « .028,  W0  » 1.6xl0*2,  (ve/c)2  «=  2 . 5xl0~* 

k*kQ  = cos  (45°).  For  this  smaller  kQ,  the  disjoint  struc- 
tures of  Figure  two  have  coalesced.  We  can  write  the  condi- 
tion for  this  coalescence  in  terms  of  P and  kQ: 

3km-  (?.  ^ 2T, 

when  - Slk«,)=0,  from  considerations  of  the  width  of 

the  unstable  regions  and  the  distance  between  the  regions  of 
instability.  Notice  that  when  , that  this  is  auto- 

matically guaranteed  to  be  the  case  by  (20)  . The  broader 
region  of  instability  is  due  to  the  larger  value  of  if  L used 
in  this  graph. 


Figure  five.  The  growth  rate  as  a function  of  k 
for  coupled  excited  daughter  waves  of  the  2«pe  instability 
driven  by  a standing-wave  pump.  This  is  a polar,  contour 
plot  of  the  solutions  to  (28) , with  the  parameters  of  Figure 
four.  The  symmetries  are  the  same  as  in  Figure  three.  The 
units  on  the  contours  are: 

0 0 

1 2x10-3  u, 

2 4x10-3 


43 


shows  a polar  plot  of  the  positive  imaginary  part  of  this 
solution. 

Thus,  the  most  important  change  that  a standing- 
wave  pump  brings  about  is  a change  in  the  effective  field 
strength  that  drives  the  instability.  As  we  have  seen 
above,  when  each  component  of  the  pump  drives  independent 
daughter  waves,  the  effective  pump  field  amplitude  is  down 
by  one-half  from  the  travelling -wave  case,  and,  when  both 
act  simultaneously,  the  amplitude  is  down  by  a factor  of  iT. 
In  this  second  case,  the  pump  fields  seem  to  add  incoher- 
ently, with  powers  rather  than  amplitudes  adding.  The 
difference  between  cases  A and  B can  be  one  of  kQ  getting 
smaller;  as  kD  decreases,  and  we  approach  a uniform  pump, 
the  effective  field  strength  increases.  The  inequalities 
in  (20)  determine  whether  the  daughter  waves  are  "broad" 
enough  to  be  pumped  by  both  components  of  the  pump  simul- 
taneously (see  Figure  one  b) , page  26). 

Indeed,  if  we  can  solve  the  full  equation  (19)  in 
the  limit  of  kQ  = 0,  we  should  recover  the  travelling- 
wave  results.  The  limit  of  (19),  however,  is  subtle. 
Strictly  speaking,  the  instability  does  not  exist  when 
kQ  = 0,  since  all  the  susceptibilities  are  proportional 
to  kQ  when  kQ  « k.  Thus,  we  can  look  at  the  limit  kQ«  k, 
but  with  kc  still  non-zero.  There  are  at  least  two  ways 
of  taking  the  limit: 

1).  In  (11), 
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Eo  * ltd [0'a)5,’w(fc#- x — o(  siV»  (k.  Z'Cj,t  *<f)]  , 

let  <X  go  to  one,  followed  by  kQ  getting  small.  This  limit 
presents  no  problems;  when  <*  = H,  it  easily  gives  the  results 
for  a standing-wave  pump,  and,  when « = 1 or  0,  it  gives  the 
properly  normalized  amplitude  for  a travelling-wave. 

2).  If,  we  setOf  = Jj,  but  take  the  limit  kQ<<  k, 
the  inequalities  in  (20)  demand  that  more  and  more  modes  be 
coupled  together.  If  we  require  that  nkQ<<  k for  any  n, 
then  we  must  use  an  arbitrarily  large  matrix  equation  of 
the  form  of  (19).  In  the  limit,  all  the  susceptibilities 
are  the  same,  and  we  have  to  evaluate  the  following  deter- 
minant, with  the  definitions: 


mi  Or)*  w.h  (a*)  (fc-S). 

[w  - S ♦ iy  ] • 

evaluated  at  the  k for  which  6/L(k)  = 8, so  that 

€ (**>-  60,)  = “ R , 


R 

0 0 

0 

Vi 

r*~ 

~ViTo“’ 

0 

0 

0 

0 

0 

0 

o Vi 

R 
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We  have  not  been  able  to  mathematically  reduce  this 
infinite  determinant  to  the  equation  for  a 2wpc  instability 
pumped  by  a travelling-wave  with  amplitude  twice  the  ampli- 
tude of  either  component  of  the  standing-wave  (this  is  the 
limit  that  would  be  expected);  however,  we  have  been  able 
on  physical  grounds,  to  derive  a lir.it  on  the  size  of 
the  determinant,  and  within  this  limit  our  large  detcr- 
minental  condition  reproduces  the  travelling-wave  results. 
That  is,  an  m x m determinant  (m  is  the  limit  we  are  able 
to  derive) , for  large  m,  gives  the  same  results  as  the 
center  2x2  determinant  (marked  by  dotted  lines  in  (35)) 
does  with  EQ  doubled.  Physically,  this  describes  the  fact 
that  as  kQ  becomes  negligible  compared  to  k,  the  two  com- 
ponents of  the  standing-wave  become  indistinguishable, 
and  the  local  amplitude  of  the  field  goes  from: 

[_to  + Ee(-^,w„)]/2  to  Ct  (wc) 

In  a straightforward,  although  tedious,  calculation 
it  is  possible  to  derive  recursive  relations  relating  the 
determinants  of  different  dimension,  m. 


(36) 

Dm  - RDm-l  " (J£/l)1Dm-t  , 

modA  m - 0, 

(37) 

Dm  = RDm-l  Dm-!  , 

mcd4  1 > 

(38) 

Dm  ""RDm-l  + PDl)  Dm-1  ) 

mod*  m ^2, 

(39) 

Dm  '-RCU-,  , 

rY)oS4n)  * 3 

where  mod4  m is  the  modulo,  base  4,  of  m (i.e.,  the  integer 
remainder  upon  division  of  m by  4).  There  are  four  rela- 
tions because  of  the  alternating  sign  of  the  diagonal 
elements  in  (35) . 

Keeping  terms  fourth  order  in  EQ,  but  no  higher 
powers,  we  can  write  the  determinental  conditions  for  the 
two  cases  of  interest,  mod4  m “ 1 or  3 (which  have  the  same 
value  to  this  order  in  LQ) : 

(40!  Dm  = [R4  * (m-0  ;‘R'/4  - (P* tfl«]  R».4  * o • 

. . 2 

It  is  easy  to  factor  this  quadratic  in  r and  write  its 

generally  complex  solutions  in  the  form: 

(41)  Dr*  01  R*  “ £ (*/2.)  • 


For  a travelling-wave  pump  we  have,  in  this  notation: 


Dm  <*  R‘  - 4 (%/zf  , 

Our  3x3  results  are  that: 

Dm  * R‘  - 2 (X/zf  ■ 

For  mN<  9,  we  get  the  expected  result,  that  j3  approaches 
four;  at  m = 9 , ft  = 3.94.  This  is  physically  equivalent 
to  the  two  components  of  the  standing-wave  field  coalescinq, 
doubling  the  local  field  amplitude  and  thus  quadrupling 
Numerical  results  shewing  this  dependence  of  f>  on  m are 
presented  in  the  next  section.  However,  after  a gradual 
asymptotic-like  approach  to  p * 4,  fi  begins  to  increase  as 


We  must  look  for  a physical  explanation  for  the 
breakdown  of  our  matrix  formalism.  Physically,  the  problem 
seems  to  in  our  continuing  to  treat  all  the  modes  as  dis- 
tinguishable as  kQ  decreases,  no  matter  how  close  they  are 

to  each  other.  That  is,  as  k decreases,  all  the  modes  at 

o 

ki2nk0  must  coalesce  with  the  mode  at  (k,  w)  and  all  the 
modes  with  wavenumbers  of  (k±(2n+l)k0)  must  coalesce  with 
the  mode  at  (k,  w -«0) . Physically,  this  will  not  be  a 
discontinuous  process,  obviously,  collapsing  an  infinity 
of  modes  for  any  non-zero  kQ  to  only  two  when  kQ  equals 
zero.  Rejecting  that  possibility,  we  must  examine  the 
criteria  for  being  able  to  distinguish  modes  which  are 
close  to  each  other. 

Referring  back  to  Figure  one  b)  (page  26) • it  is 
clear  that  if  kc  is  small  enough  so  that  for  n»l,  we  have: 

} £0a(kj)- CA^aO?a-2h^i  « 

and,  then,  in  this  case,  it  will  be  difficult  to  distinguish 
between  t*>a  ( ka)  and  t^(kai2k0).  Both  of  these  modes  will 
overlap  entirely  and  so  must,  physically,  be  considered 
identical.  This  result  puts  a rough  upper  limit  on  n: 
n 4-  10,  and  on  m,  the  size  of  the  determinant,  m^  20.  This 
consideration  gives  an  upper  bound  on  the  number  of  modes 
which  can  be  put  in  this  matrix  formalism  of  (19)  and 
still  remain  physically  meaningful.  That  is  to  say,  for 
m large  enough,  say  greater  than  10,  we  expect  the  linear 
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results  for  these  modes  to  be  meaningless,  because  the 
modes  are  no  longer  distinguishable.  And,  it  is  at  this 
point  that  our  mathematical  results  begin  to  give  physically 
suspicious  results.  Another  way  of  looking  at  this  situa-' 
tion  is  to  say  that  only  on  the  order  of  ten  modes  with  the 
same  frequency  can  be  excited,  and  the  coupling  of  these 
ten  modes  is  given  by  (34).  When  this  many  modes  must  be 
taken  into  account,  we  have  the  standinq-wave  pump  duplica- 
ting the  coupling  given  by  a travelling-wave  pump. 

It  is  important  to  note  that  this  argument  is  a 
physical  and  intuitive  one.  Mathematically , in  an  infinite, 
homogeneous  plasma,  there  are  no  limits  on  the  number  of 
nodes  that  can  exist.  If  tv/o  modes  have  different  wavenum- 
bers, no  matter  how  close  these  wavenumbers  may  be  to  each- 
other,  they  are  distinguishable.  For  a reason  that  remains 
unclear,  our  matrix  formalism  does  not  mathematically 
allow  us  to  take  the  limit  kQ-*0  and,  by  doing  so,  to 
reproduce  the  travel ling -wave  results. 

THE  ELECTRON-ION  DECAY  INSTABILITY1'2 

In  the  case  of  the  electron-ion  decay  instability, 
the  decay  of  a high  frequency  electric  field  near  the  plasma 
frequency  into  a low  frequency  acoustic  mode  (or  quasi-medo) 
and  a high  frequency  Langmuir  wave,  the  new  mode  coupling 
effects  of  a standing-wave  pump  are  similar  to  those  for 
the  2WpC  instability.  The  important  differences  between 
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the  standing-wave  treatment  of  these  two  instabilities 
come  from  the  fact  that  the  daughter  modes  for  the  decay 
instability  are  different  resonances  of  the  system.  Thus, 
to  begin  with,  we  must  interpret  the  inequalities  in  (20) 
differently.  The  number  of  modes  of  the  acoustic  branch 
which  are  close  enough  (in  the  sense  above)  is  na: 

(43  ) 2.  Cry  0 (Wt/w,/1  k,  < (uy 

and  of  Langmuir  waves,  is  nj-,: 

<44>  t ^ kp  '(k  ♦hbke)  > 6 fa,  k ♦AfeJ  • 

In  the  case  where: 

(45)  -3a*  +3kcc$&  > (r'u/m.)*'*,  & = Ccr,~'(k- k*) , 

these  relations  imply  that  the  standing  wave  pump  will 
couple  together  first  modes  at: 

c Oi(k ),  Ua(4  + &o),  (/)A(k-k„), 

and  then  additional  modes  as  allowed.  When  (45)  holds, 
there  cannot  be  a coupling  of  one  acoustic  wave  and  two 
Langmuir  waves.  This  follows  from  the  frequency  matching 
condition: 

<Vo  ~ CVa  + , 

which,  for  a transverse  pump,  assures  that  k is  always 
much  greater  than  k . Thus  for  any  value  of  k_,  n3  will 
be  greater  than  nb. 

We  have  looked  in  some  detail  at  this  case,  the 
simplest  standing-wave  coupling,  that  of  wL(k)  coupled 
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with  ‘-'A(k±k0)  = Wo-(0L(k).  Corresponding  to  equation  (28) 
for  the  2«pe  instability,  we  have  (from  (25)): 

tu-  w-  (k  '■k.)  * i%i][  Lu-  oj,  * U)*  Ck-k.)  * 

(46)  -<xz  (V 

-(.i-*)1  rl  [ to- * o , 

where  (see  (A. 8)  and  (A. 9)): 

(47)  iV  ' ' -Zf  (^)^o  (&&)z  ■ 

This  equation  is  valid  as  long  as  (using  (43)  and  (44)), 
and  na  = n^  = 1 : 


(48a) 


CWc/W;)^  k0  < 


(48b) 


V*n  € (c o)  < 6A-A#  + 6 k»z. 


Notice  that  as  long  as  9t*  0. , then  = 0 (me/ma)  ljk)  , so 
that  (48a)  is  easily  satisfied;  (48b) , however,  is  a re- 
strictive condition  on  Im€(WL).  It  is  always  the  case, 
furthermore,  that: 

(*<>.) | « (*>o  . 

This  means  there  is  no  regime  corresponding  to  case  A for 
the  2wpe  instability,  and  it  is  always  a valid  approximation 
to  write  for  (4  6),  when  ot* 

(49a)  - Z(r/Z)1-Q. 
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That  is,  for  the  decay  instability  with  a standing-wave  pump 


and  kQ  in  the  range  prescribed  by  (48a)  and  (48b) , the 
coupling  is  half  what  it  would  be  for  a travelling-wave 
pump  with  the  same  local  amplitude.  The  effective  power 
is  down  by  a factor  of  two.  The  formal  solutions  to  this 
type  of  equation  (49a)  have  been  studied  in  detail  else- 
where, and  the  growth  rate  is  gotten  straightforwardly 
by  solving  (49a) : 


whe  re : 


A = ^A(k)  *•  GJl  (k)  - Uj0 


This  oives,  with  4=  0: 


(49b) 


IV 


max 


«*a  , &«  yA  , 

- - t r’»*A  , Pl»V.  , 

. I"/**  , *Uj  *,.»#*. 


We  note  here  only  the  consequences  of  a standing-wave  pump. 
These  changes  all  result  from  the  factor  of  two  in  (49a). 
This  halving  of  the  effective  power  means  that,  given  a 
local  field  amplitude,  the  threshold  is  doubled  if  this 
field  is  a stand ing-wave  rather  than  a travelling  wave.  In 
the  case  of  a travel ling -wave  with  local  amplitude  LQ,  the 
threshold  for  the  decay  instability  is: 


(50a) 


(Tc/r)1  = k a. 


f 
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and  the  maximum  growth  rate  is: 


u>  | - 

\ , 

(50b) 

r/i  , 

r » 

. 0% , 

r«*L 

, & 

But  if  the  pump  field  is  a standing-wave  with  total,  local 
amplitude,  r.Q,  then  the  threshold  is  given  by: 


(Cc/i)1  = 2.  u* 


and  the  maximum  growth  rate  is  given  by  (49b) . 

Let  us  look  at  one  situation  where  this  factor  of 
two  may  be  important.  In  several  ionospheric  modification 
experiments,  it  has  been  claimed  that  the  observed  enhanced 
level  of  Langmuir  waves  is  due  to  a non-equilibrium  spectrum 
created  by  the  decay  instability . ^ ' *0-15  Physically,  the 
scenario  is  as  follows:  the  relatively  intense  electro- 

magnetic wave  from  the  ground  (a  wave  of  ordinary  polariz- 
ation) travels  upwards  from  the  antenna  through  the  iono- 
sphere. The  inhomogeneity  of  the  ionosphere  has  three 
effects:  it  results  in  an  Airy  swelling  factor  in  the  ampli- 

tude of  the  modifier  wave,  in  an  absorption  which  becomes 
a function  of  altitude,  and  in  a reflection  of  the  modifier 
wave  when  the  critical  density  is  reached  (when  Wpe (height) 

3 WQ) . Since  this  absorption  is  usually  small,  the  modifier 
wave  has  sufficient  power  to  exceed  the  decay  instability 
threshold  at  the  reflection  point  and  at  altitudes  somewhat 
lower.  If  this  threshold  is  exceeded,  part  of  the  energy 
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from  the  modifier  wave  is  transferred  to  Langmuir  turbulence 
by  the  E1D  instability.  The  instability  is  driven  by  the 
sum  of  the  incident  and  reflected  waves,  a standing -wave. 

The  saturation  of  this  instability  then  provides  a steady 
state  level  of  Langmuir  turbulence  that  is  measured  from 
the  earth  using  a second  wave  (which  is  scattered  by  this 
enhanced  turbulence  in  the  "modified"  ionosphere). 

Attempts  to  understand  the  experimental  data  on 
this  scattering  revolve  around  estimating  the  ratio  of 
modifier  intensity  to  the  LID  instability  threshold,  then 
calculating  the  enhancement  of  the  Langmuir  spectrum  that 
would  be  created  by  such  a pump,  and  finally,  calculating 
the  scattering  that  such  an  enhanced  spectrum  would  pro- 
duce. Our  results  will  fc--  relevant,  since  the  modifier 
wave,  is,  in  fact,  a standing-wave  whenever  it  is  of 
sufficient  intensity  and  satisfies  the  frequency  matching 
condition  for  the  instability.-*  in  a detailed  study  of  the 
problem,  DuBois,  et  al.  show  that  even  with  considerable 
absorption  of  the  incident  wave,  a standing-wave  field 
is  set  up  considerable  distances  from  the  reflection  point. 
Using  the  parameters  of  this  paper,-*  which  correspond  to 
a height  about  560  meters  below  the  reflection  point,  we 
have:  kQ  = 4x10”’’  (in  units  of  the  Debye  wavenumber,  which 

is  here  220m-*-),  the  fastest  growing  mode's  wavenumber, 
k = 0.05,  and  k*£0  = 0.85.  In  addition,  DuBois,  et  al,, 


calculate  a local,  average  power  flux  of  800  W/m^  at  this 
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height,  for  an  antenna  inensity  flux  of  50  W/m2.  This  is 
about  five  times  the  threshold  of  the  EID  instability,  so 
that : 

Wv  € M — 2.  (jj  ta  2»|0'r 

(in  units  of  the  plasma  frequency).  Thus,  inequality  (48b) 
is  satisfied.  Notice,  however,  that  even  for  a slightly 
higher  power,  (48b)  would  not  be  satisfied,  and  we  would 
have  to  solve  a more  complicated  coupling  of  waves. 

Thus  equation  (49)  describes  the  decay  instability 

in  this  case.  The  correct  way  to  account  for  the  reflected 

wave  in  the  ionosphere  is  not  to  add  amplitudes,  but  rather 

• 

to  add  fluxes  or  powers  of  the  two  components  of  the  standing 
wave  (that  is,  consistent  with  our  homogeneous  approximation, 
we  must  "decompose"  the  actual  Airy  function  amplitude 
modulation  of  the  modifier  wave  and  look  at  it  as  made  up 
of  two  waves,  with  different  relations  betweeen<^  and  kQ 
for  each).  Thus,  the  estimates  of  the  effective  power  in 
the  papers  referred  to  above  are  too  high  by  a factor  of 
two.  Even  given  the  large  uncertainty  in  all  the  calcul- 
ations of  the  ionospheric  problem,  a increase  by  a factor 
of  two  in  the  threshold  is  critical.  Perkins,  et  al.,  for 
example,  estimate  the  field  at  a value  only  twice  the 
threshold, *4  and  DuBois  and  Goldman  use  a figure  of  five 
times  threshold. 3 in  the  graphs  which  DuBois  and  Goldman1-5 
show  of  total  (saturated)  Langmuir  energy  as  a function  of 
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pump  power,  a change  by  a factor  of  two  in  the  effective 
pump  power,  from  5.0  to  2.5,  results  in  a decrease  in  the 
peak  Langmuir  wave  energy  by  a factor  of  100,  and  in  the 
integrated  Langmuir  energy  by  a factor  of  20.  Since  all 
the  calculations  of  the  scattering  that  should  be  observed 
are  several  orders  of  magnitude  too  large,  this  changed 
threshold  helps  in  accounting  for  the  discrepancy.  Cer- 
tainly, this  is  not  a conclusive  demonstration  of  the 
effect  of  a standing-wave  pump,  but  our  interpretation  of 
the  ionospheric  modification  experiments  does  open  up  a 
means  of  testing  the  effect  of  a standing  wave  pump.  This 
factor  of  two  may  be  observable  in  the  comparison  of  low 
power  modification  experiments  (with  flux  densities  of  less 
than  50  W/m^) , which  have  the  simplest  coupling  of  three 
daughter  waves  described  above,  and  higher  power  experiments, 
in  which  (48b)  is  no  longer  satisfied  and  more  modes  are 
coupled.  The  theory  of  the  case  when  kQ  is  smaller,  or  the 
power  greater,  so  that  (48b)  is  no  longer  satisfied,  is 
treated  next. 

As  in  the  case  of  the  2A>pe  instability,  we  have  also 
looked  at  the  problem  of  the  travel ling -wave  limit  of  the 
decay  instability.  Rather  than  an  analytic  treatment,  we 
have  numerically  solved  the  equivalent  of  equation  (34)  for 
the  decay  instability.  Specifically,  we  tested  each  mode, 
at  the  k for  perfect  matching,  to  determine  whether  it  was 
"close  enough"  to  u\(k)  to  be  coupled.  Thus,  the  results 
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of  our  calculations  were  based  on  Table  III,  where  the  * 

number  of  modes  of  the  form  wv(k  t nkQ)  , b0  - Oa  (k  t nkQ) 
that  are  within  Anw  of  LMk)  are  listed.  Thus,  coupling 
three  Langmuir  modes  means  that  there  are  26  acoustic  modes 
that  are  close  enough  to  be  coupled.  We  actually  evaluated 
a 29  x 29  matrix  equation  in  this  case,  rather  than  taking 
advantage  of  the  inequalities  (43)  and  (44)  which  would 
have  allowed  us  to  neglect  22  of  these  acoustic  modes. 

This  test  cf  the  assumptions  in  (43)  and  (44)  showed,  in 
fact,  that  these  are  good  assumptions.  Figure  six 
shows  the  results  of  this  numerical  work,  where  we  have 
graphed  the  growth  rate  of  the  instability  as  a function 
kD.  Here  the  growth  rate  for  k0  = .006  has  the  value  from 
by  equation  (49).  As  more  modes  are  coupled  together,  the 
growth  rate  increases,  it  approaches  the  value  that  would 
obtain  for  the  solution  to: 

(co-d4a).ijf)(to-co. - r*  =o, 

exactly  the  result  for  a travelling  wave-pump  with  the 
same  local  amplitude  as  the  sum  of  the  amplitudes  of  the 
two  components  of  the  standing-wave  pump.  Figure  seven 
presents  the  same  results  as  a function  of  the  number  of 
modes  coupled  together.  This  shows  the  dependence  on  the 
smaller  of  na  and  that  we  noted  above;  it  is  a result 
of  our  matrix  equation  which  allows  no  coupling  between 
waves  which  have  wavenumber  shifts  which  differ  by  more  than 


TABLE  I II 


MODES  COUPLED  IN  THE  DECAY 
INSTABILITY  AS  ko-*-0. 


ko 

2na, 

Number  of 
Acoustic  Modes 

2nb-l, 
Number  of 
Langmuir  Modes 

^max 

.0062 

1 

1 

.76x10"? 

.0031 

2 

1 

1 . 20x10*4 

.00041 

8 

1 

1 . 34x10*  ^ 

.00020 

14 

1 

1.33x10*4 

.00016 

18 

1 

1.34xl0“4 

.00014 

22 

3 

1.79x10-4 

.00012 

26 

3 

1.81x10-4 

NOTE:  na  and  nb  and  the  number  of  modes  are 

calculated  according  to  equation  (19).  The  limit  as  kj*0 
was  taken  keeping  0>o  constant,  in  this  case,TL"  10~4 
and  £ « . 6 . 
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Figure  six.  The  growth  rate  of  the  electron-ion 
decay  instability,  pumped  by  a standing-wave  electro- 
magnetic waver  as  a function  of  the  pump  vavenumber. 

For  this  figure,  WQ  = 4.xl0“7,  (me/mi)  ■ 5.xl0”3.  The 
growth  rate  increases  as  the  number  of  modes  coupled  to 
gether  increases. 


Figure  seven.  The  growth  rate  of  the  electron-ion 
decay  instability,  pumped  by  a standing-wave  electromagnetic 
pump,  as  a function  of  the  number  of  waves  that  are  coupled 
together;  the  parameters  here  are  the  same  as  in  Figure 
six. 
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2kQ.  In  Figure  eight,  we  have  fitted  the  following  equation 
to  the  growth  rates: 

(52)  £fc>-Wt( fiMy)  (to + 0- 

and  graphed  ft  against  the  number  of  modes.  The  analytic 
description  of  this  limit,  derived  for  the  2fc»pe  instability, 
is  also  applicable  here.  Note,  for  example,  that  in  the 
case  that  m = 7,  here  corresponding  to  three  Langmuir  waves 
coupled  with  four  acoustic  waves,  equation  (35)  predicts 
that  ft = 3,  as  our  numerical  work  confirms. 

Finally,  we  note  that  when  k*kQ  = our  analysis 
(see  equation  (44))  predicts  that  the  standing-wave  effects 
will  be  negligible  whenever,  in  addition: 

(53)  6 K»  -£>71  € (cc ) 

Physically,  this  is  a consequence  of  the  fact  that,  given 
k»  kQ,  the  variation  of  Wx,(kikp)  with  kD  comes  almost 
entirely  from  the  k*kD  term: 

3&  /L.J  ■ 

If  k«kc  = 0,  all  the  Langmuir  waves  at  ^(kinkjo)  will  be 
very  close  together,  and  ft-* 4 as  these  modes  are  coupled. 


Figure  eight.  "Effective  pump  power,"  ft  , as  a function 
of  the  number  of  modes  coupled  together,  ft  is  defined 
in  equation  (52).  4 corresponds  to  the  amplitude  of 

the  field  contributing  as  in  the  case  of  a travelling-wave. 
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CHAPTER  III 


LANGMUIR  PUMPED  PARAMETRIC  INSTABILITIES  IN 
THREE  DIMENSIONS 


In  this  chapter  we  will  study  the  four-wave  pro- 
cesses which  have  longitudinally  polarized  Langmuir-like 
waves  as  Loth  pump  and  daughter  waves.  These  instabili- 
ties are  well-known  in  one-dimension  (for  references  sec 
Table  I,  page  6,  and  II,  page  8),  when  the  daughter 
waves  have  wavevectors  parallel  to  k^j.  These  well  known 
one  dimensional  instabilities  for  the  Langmuir-to-Langmuir 
case  are  tne  oscillating  two  stream  instability  (OTS) , 
which  couples  the  high-frequency  daughter  waves  with  a 
zero  frequency  response  of  the  plasma,  and  the  modified 
EID  instability,  which  couples  the  high  frequency  daughter 
waves  with  a resonant  ion  acoustic  mode.  In  our  analysis 
of  these  instabilities,  we  have  derived  a single  dispersion 
relation  wich  includes  all  these  instabilities,  and  which 
includes  regimes  where  these  instabilities  exist  for  wave- 
numbers  other  than  those  parallel  to  kQ.  We  have,  in  ad- 
dition, discovered  another  instability,  the  stimulated 
modulational  instability  (SM)  which  only  appears  in  three 
dimensions.  Our  treatment  for  these  instabilities  in- 
cludes the  effects  of  an  arbitrary  wavevector  pump  and 
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an  exact  ion  equation  (rather  than  the  usual  fuild  treat- 
ment) . 

As  the  kinetic  theory  derivation  of  the  last 
chapter  implied,  we  can  express  the  susceptibility  of  the 
plasma  in  a power  series  in  the  fields,  so  that,  at  least 
formally,  we  have: 

a)  x(e)  = &-F  E r£  (£  45  §)+  " 

Physically,  based  on  an  analogy  with  quantum  mechanics, 
the  higher  order  terms  correspond  to  the  plasma  couplina 
with  more  than  one  pump  field  photon.  In  the  last  chapter, 
by  restricting  our  work  to  phenomena  involving  only  we 
were  able  to  deal  with  only  two  daughter  fields,  a 
plasmon  and  photon,  two  plasmons,  or  two  photons.  However, 
important  new  effects  occur  with  the  couplings  involved 
in  ^3,  the  so-called  four-wave  interactions.  To  treat  the 
these  interactions,  we  must  include  effects  due  to  both 
and  arbitrary  kQ,  since,  in  the  case  of  a Langmuir  or 
Langmuir-like  wave  (for  both  pump  and  daughter  wave) 
one  has  a dispersion  relation  which  allows  daughter  waves 
with  wavenumbers  comparable  to  the  pump  wavenumber.  This 
follows  since  the  dispersion  of  a Langmuir  wave  is  pro- 
portional to  ve2  rather  than  c2,  as  in  the  case  of  an 
electromagnetic  wave.  The  larger  pump  wavenumbers  make 
three  dimensional  effects  important.  After  a theoretical 
study  of  the  linear  stability  of  a Langmuir  pump  to 


our  results  to  the  Langmuir  spectra  associated  with  electron 
beams  in  several  parameter  regimes:  the  solar  corona,  the 

ionosphere,  and  low  density  laboratory  plasmas.  Our 
general  conclusion  from  these  applications  is  that  para- 
metric processes  do  not  greatly  affect  the  dynamics  of 
the  beam  plasm  interaction  (as  has  been  theorized) , but, 
even  so,  that  these  instabilities  can  account  for  a number 
of  other  phenomena  associated  with  these  beam  plasma  systems. 

DERIVATION  OF  THE  DISPERSION  RELATION 1 

Using  the  formalism  developed  thus  far,  it  is 
straightforward  to  write  down  the  determinental  cone  .tion 
describing  the  parametric  processes  induced  by  a pump  made 


up  of  Langmuir  waves,  that  is,  a pump  that  is  high  fre- 
quency and  longitudinally  polarized,  with  frequency  and 
wavenumber  k^,  (w  is  a low  frequency  for  all  of  this  chapter)  : 


f € (H,tv) 

/e(K,w)  ^ 

, 

(1) 

€(k-kf/io-cJ 

.Vwi;  (?•'£*, in  Ec 

E(k-k», 

1 

= 0, 

6 / 

EtVki,  , 

f 

/ 

^♦h)  1 

where,  the  entries  in  the  first  row  and  first  column  are 
the  susceptibilities  derived  in  the  Appendix  from  the 
second  order  current,  and  the  entries  in  the  (2,3)  and 


(3,2)  positions  are  the  susceptibilities  derived  from  the 
thrid-order  current.  These  higher-order  susceptibilities 
also  enter  into  the  definitions  of  the  diagonal  elements, 
"renormalizing"  the  dielectric  functions: 

(2)  e(M  - £(*,<«>)  4 /-(&*>;•!!•  & ' 

Notice  that  the  upper  2x2  matrix  gives  the  usual  decay 
instability,  in  which  case,  as  is  well-known,  the  higher- 
order  susceptibilities  are  not  necessary.  However,  when 
both  the  Stokes  and  anti-Stokes  modes  are  coupled  in  the 
problem,  they  are  each  coupled  (in  first  order)  to  the 
low  frequency  field  at  E(k,u>)  and  (in  second  order)  to 
each  other.  Since  we  are  dealing  with  Langmuir  oscil- 
lations, these  matrix  elements  must  be  evaluated  for  ar- 
bitrary ko?  we  cannot  use  the  previously  published  results 

, A A 

which  assumed  either  k >>  kQ  or  k*kD  - 1. 

The  rather  tiresome  derivation  of  the  susceptibil- 
ities is  deferred  to  the  Appendix.  Even  without  the  ex- 
plicit form  of  these  elements,  we  can  use  the  following 
relations  to  write  the  determinate  of  (1)  in  a strikinqly 
simple  form  (see  (A. 8),  (A. 10),  and  (A. 11)  from  Appendix): 

(3)  k,, I*,, fet uttt) r fe  (^t ^ , 

and : 

^ ( k*  *k.  tw,i  m k1  l/(  k*fc,  Wtw,;  t£,,*  |*j  K,  w)  • 

/(*%***;  *&'.**+,■  &*)  I • 


(4) 


With  these  relations,  all  terms  third-order  or  higher 
in  Eq  cancel  and  we  get,  after  rearrangement: 

I k'e  Ck,  X(k‘h.  u-tu>c;  >(’•■/, , uJtUc)  E-J 

' kx  (.(*,(-*)  -I  J 6(k+£,,UtV>) 

(5) 

t XflNr.. 

£Qi-kt'U>-De)  o 

This  equation  shows  very  clearly  the  coupling  of  the  three 
resonances  of  the  plasma,  two  with  high  frequency  and  one 
with  low  frequency,  in  a form  characteristic  of  mode- 
coupling in  non-linear  optics,  for  example.  Using  the 
explicit  expressions  for  the  susceptibilities  (equations 
A. 8,  A. 10,  and  A. 11,  from  the  Appendix),  we  have: 


t J kz£(!s,v)  ) w.  r MS 

/ 4 Lec *♦*(.,  ww0) 

6(k-ks^-^c)  j 

where : 

T*  0*/©iv 

r rk  EoMfe.) -E.(fe)  . 

(7)  ne  (0C  t &;)] 

A'5  [k-(k*fe)]2/[fc  |J?.i  kjT  ■ 

Several  aspects  of  equation  (6)  should  be  noted: 

1) .  The  two  high  frequency  waves,  the  daughter 
(Langmuir- like)  waves  at  k t J^,  are  coupled  together 
through  a low  frequency  response  of  the  plasma.  As  Tsyto- 
vich  noted  in  his  book,  these  couplings  dominate  all  other 
Langmuir  wave  interactions.2 

2)  . This  low  frequency  response  of  the  plasma  is 
not  given  by  the  plasma  dielectric  function,  except  in  the 
limit  that  f-(k,  w)  goes  to  zero.  Let  us  call: 
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(3) 


D(fc,w) 


T kX€  (k  tQj) 

T>i"  kz €^“»  ■ 


At  frequencies  for  which  k2 €(k,  w)  is  not  sm?il  comapred  to 
one  (k  is  in  units  of  k ),  the  plasma's  low  frequency 
response  is  complicated.  This  "compensation"  has  been  noted 
before,-*'4  and  arises  from  the  terms  that  the  third-order 
current  contributes.  Figure  one  shows  "diagrams"  of  the 
contributions  to  the  wave-wave  interactions  from  the  second- 
order  current  (in  Figure  one  a))  and  the  dominant  contri- 
bution from  the  third  order  current  (in  Figure  one  b)). 

The  terms  that  contribute  from  the  third  order  current 
are  essentially  second  order  effects  coupled  together  by 
a low  frequency  response  of  the  plasma;  the  "four-wave" 
diagram  can  be  viewed  as  two  of  the  three-wave  diagrams 
joined  by  a (virtual)  low  frequency  response. 

3) . Equation  (6)  also  comes  out  of  a more  careful 
analysis  of  the  linear  stability  of  an  arbitrary  Langmuir 
spectrum,  if  that  spectrum,  centered  about  u,0(k0),  has  a 
spread,  &k0,  in  k-space  about  kQ,  such  that  a condition  on 
the  narrowness  of  the  spectrum  is  satisfied: 


(9) 
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This  condition  insures  that  the  resonance  is  a broader 
function  in  k-space  than  the  spectrum,  so  that  the  detailed 
shape  of  tne  spectrum  can  be  ignored. 


Fiaure  one*  Three-wave  (a)  and  four-wave  (b) 
interactions.  The  pump  wave  has  wavevector  and 
frequency  DQ;  the  low-frequency  electrostatic  disturbance 
has  wavevector  k and  frequency  w;  and,  the  Langmuir 
daughter  waves  have  wavevectors  kQ  i k and  frequencies 
b0  ± w.  In  the  four-wave  case  both  oF  these  hiqh 
frequency  waves  are  excited  and  in  the  three-wave  case, 
only  one  of  them, the  Stokes  wave,  at  k_  - k,  Q0  - u . 

The  wave  at  Cjq  + w is  called  the  anti^tokes  wave. 
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The  question  of  the  linear  stability  of  a Langmuir 
pump  has  been  treated,  in  various  approximations  (although 
never  before  in  three  dimensions),  by  many  authors;1,4-8 
our  results  generalize  these  previous  treatments  in  throe 
respects:  we  have  ^ arbitrary;  we  have  treated  the 

problem  in  three  dimensions,  with  cylindrical  symmetry 
about  k^;  and,  we  have  used  an  exact,  kinetic,  expression 
for 

do)  k'eflj.w)  — I +r  + (k^a)Z(j^ff)  + T Z (k^f) , 

where  Z is  the  Fried-Conti  plasma  dispersion  function.9 
We  can  connect  our  results  with  previous  theoretical  treat- 
ments which  used  the  following  "fluid"  form  for  D(k,  u>)  : 

ui)  W + , 

where  is  a phenomenologocal  damping  of  an  ion-acoustic 
mode  or  quasi-mode  at  u=  Ufa- iVA/2.  Our  exact  equation  for 
D(£.,  w)  reduces  to  (11)  when: 

ftV;  4<.  60  kv«  , 

and  we  then  can  use  the  asymptotic  forms  for  the  Z-function 

\ « 1 ; 

I* ~ip "4*  ’ ^ ■ 

Substitution  of  these  forms  for  Z into  (10)  and  the  defini- 
tion of  D (k , cj)  gives  the  real  part  of  (11);  the  imaginary 
part  can  also  be  gotten  similarly.  However,  the  fluid 


f 
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treatment  qoes  on  to  use  (11)  for  all  values  of  W.  A speci- 
fic problem  arises  with  the  zero  frequency  limit  of  (11). 

For  v=  0,  (11)  gives  the  value  -1.  This  is  also  the  value 

of  the  static  limit  of  the  exact  expression  for  D(k,  w)  . 

In  fact,  (11)  and  the  exact  form  for  D(k,  w)  agree  only 
when  = 0 and  w=  0(wA).  As  we  shall  see,  this  is  not  a 
severe  limitation  on  the  fluid  expression,  since  it  is 
only  in  these  two  regimes  that  interesting  solutions  exist 
for  our  dispersion  relation,  (6) , and  quantitatively  satis- 
factory results  can  be  obtained  with  a suitable  choice  of 
WA  and  vA. 

In  the  calculations  presented  here,  we  have  made  a 
resonance  approximation  for  the  high-frequency  dielectric 
function: 


1 


J 


d2)  eCkt^vtub)  - ±2-[oj±uc  TvL(k±k*)*  iyL]  , 

being  phenomenological  damping  (either  collisional  or 
Landau)  and  an  exact  expression  for  the  low  frequency 
C(k,  Co)  in  terms  of  the  Fried-Conti  Z-function  given  in 
equation  (9).  Finally,  in  the  frequently  valid  approxi- 
mation that: 


(which  is  usually  guaranteed  because  **  yu_)  , (6)  becomes: 

(i3)  D(k,w)[(W“X#  MyJ2-"Sl]+  (/*ss.  vCS*)  -o, 


where : 


= (ktko)  = TX0  , 

25  ^ Zuo-uXk^-UiCk-k',)*-  3k1, 

2/o  - 5_ ' S+  - 6 k- fee  , 

in  dimensionless  units  of  wavelength  in  units  of  kD  -1  and 
frequency  in  units  of  the  plasma  frequency.  is  the 

cosine  of  the  angle  between  the  pump  wavevector,  J^,  and  the 
wavevector  of  the  duaghter  waves,  kik^,.  The  three  deltas 
are  different  frequency  mismatches:  8+  being  the  mismatch 
between  the  up-  (or  down-)  shifted  wave  and  the  pump;  and 
S being  the  average  of  these  mismatches. 

Before  proceeding  to  applications  of  equation  (6) , 
we  need  sane  way  of  systematizing  the  large  number  of 
regimes  in  kQ,  WQ,  VA,  and  in  which  qualitatively  differ- 
ent solutions  to  (6)  exist.  We  will  restrict  our  results 

i 

in  tnis  section  on  analytic  results  to  the  case  that 
This  is  almost  always  satisfied.  We  have  organized  our 
study  of  the  analytic  solutions  to  (6)  as  follows:  we  first 
divide  the  solutions  with  positive  imaginary  parts  into 
three  groups  on  the  basis  of  the  physical  processes 
which  account  for  the  unstable  solutions.  Each  of  these 
physical  processes  has  a group  of  parameters  which  deter- 
mine its  qualitative  features. 

1).  Oscillating  two-stream  instabilities  (OTS). 

There  are  a class  of  growing  solutions  to  (6)  character- 
ized by  their  non-oscillatory  nature,  Rew<In  u.  By 


I 

( 


analogy  with  the  well-known  OTS  instability , 4 * 10  we  call 
these  OTS-like  instabilities,  even  though  we  will  include 
a much  larger  class  of  behavior  under  this  name.  This  is 
a four-wave  interaction,  schematically  shown  in  Figure  one. 

2) .  Electron-ion  decay  instabilities  (EID) . These 
solutions  occur  when  U)^  uA  , and  only  one  of  the  high 
frequency  dielectric  functions  is  near  resonance.  This  is 
a three-wave  coupling,  and  is  shown  in  Figure  one.11*12 

3) .  Stimulated  modulational  instabilities  (SM ) . Also 
known  as  a "modified  electron-ion  decay  instability,"  this 
is  a four -wave  interaction,  in  which  both  high  frequency 
waves  participate.  Several  special  cases  of  the  SM  insta- 
bility have  been  studied  for  an  electromagnetic  pump.6'13 

OSCILLATING  TWO  STREAM  INSTABILITIES 

The  character  of  this  class  of  instabilities  is 
sufficiently  determined  by  WQ  and  kQ.  For  the  OTS  insta- 
bility, both  k and  Ira  to  depend  on  WQ,  and  we  have  been 
able  to  solve  (6)  analytically  in  the  following  regimes: 

A.  |co|  4<  U)t  , or  Vlo  «.  ; W c kfi 

0.  M «•  K , k »fc.;  or  We  «64vr  , Wc  >;  10  k cz; 

c.  M » ty  , k or  W0»64’/r  t W.  <<•  lOfc1; 

D.  M » v,  . Or  W0  »64Vr  , » I CK*  ; 


C(E  y CWe/rwi)CT’t  'JjT 

0J-,  = (»•/»*, )*  Vy 


whe  re : 
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In  the  intermediate  regions,  where  k is  on  the  order  of  kc 

- V 

and  |u»|  on  the  order  of  td/  * numerical  solutions  are  found 

which  connect  the  analytic  solutions. 

P)  j | « j ) k <<  k,  . Under  these  conditions: 

£Us,v)  * £ (K,o)  = H -~*—r  ' 

Kp,  * 


D CIs,  o)  — - ( , 


so  that  (6)  becomes: 


d4)  Cw +7.-X.)z-S2  +A,§*)  "G  ' 

If  we  call  w = x + iy,  then  an  OTS  solution  exists  when 


k-kQ  «•  k/kQ,  and 


(15)  ^ = • 

The  threshold,  maximum  growth  rate  and  wavenumber  for 
maximum  growth  are  easily  gotten;  solving  for  y = 0 in 
equation  (15)  gives  the  threshold  value  for  VJ0 

W.‘  - s/ • 

To  get  the  maximum  growth  rate  by  differentiation  of  (15) , 
we  need  to  deal  with  the  angular  dependence  in  ft*  We  can 
write  the  angular  factor  in  (15)  as: 

/tfS- ♦/*-*««  _ , ^ -•«■>  ffc1/  J + W+W* 

/*  ~ " 5 ^ [ k,1  It  + +WtS- 4k*c«vv  k.r / ) ’ 


} * C05w 


where : 


In  the  limit  that  k^<<  kQ^ : 

/<*  - I -k'sin'p/kJ 

and  differentiation  of  (15)  with  respect  to  X gives  a con- 
dition for  maximum  growth  independent  of  k:  /-~90* 

Now,  differentiating  (15)  with  respect  to  k gives  the  k for 
which  maximum  growth  occurs: 

U6a,  S.-} 

and 

U6L,  J-H-  ” 

In  fact,  differentiation  of  the  general  angular  dependence 
of  /t  shows  that  the  growth  rate,  in  all  regimes,  has 
extrema  at  f-  = 0°  and  90°.  This  is  a general  result  as  long 
as  the  growth  rate  has  its  only  explicit  angular  dependence 
through  /<.  Our  approximation  requiring  that: 

J'tn  (aj  > Rjt,  <<  OJ; 

demands  that  we  choose  the  maximam  at  90°.  The  original 
assumption  that: 

M « U){ 

gives  the  strictest  upper  bound  on  W0  and  the  lowest  bound 
on  kQ;  for  this  case: 

|5f/l  « xVl  , 
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or: 

« 64-oi/f  , 

so  that: 

1 0 I:*  «'  W0  ^ • 

Figure  two  a shows  numerical  solutions  of  (6)  in  this 
regime.  The  OTS  instability  is  the  small,  roughly  cir- 
cular feature  perpendicular  to  the  axis  at  the  point 

*o-*“llo- 

B)  . 1 | <<  0Jf  f h >>  . In  this  case,  equation  (14) 

also  describes  our  results;  however,  when  k » kQ,  then  we 

A A 

must  have  k*kQ  ss  1 (this  is  the  other  maximum  of  yu  noted 
above).  Thus,  while  x and  y are  again  given  by  (15),  the 
maximum  growth  occurs  on  axis.  For  this  regime  to  be  ap- 
plicable : 

for  the  OTS  generally.  This,  combined  with  the  other 
conditions  gives: 

W,  6+*/r  t |0  k,1  « W,  . 

Figure  three  shows  a case  not  strictly  in  this  para- 
meter regime,  but  which  has  the  same  qualitative  features. 
The  k for  maximum  growth  is  given  by  (16a) r and  we  have 
the  figure  eight  structure  of  the  previous  regime  now  on 
axis  rather  than  perpendicular  to  the  axis.  More  inter- 
esting is  the  transition  from  the  k » kQ  case,  where 
growth  occurs  for  k*k0  - 1,  to  the  k <<  kQ  case,  where  it 
occurs  fsr  £*k0  ■ 


0 
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Figure  two  a) . Solutions  to  the  dispersion  re- 
lation, equation  (6),  showing  contours  of  constant  growth 
rate,  Im  <j,  as  a function  of  the  wavevector,  k^-k,  of  one 
of  the  excited  daughter  Langmuir  waves.  Here,  WQ  = 4xl0-6. 
The  pump  wavevector  is  shown  with  an  arrow,  with  magnitude, 
kQ  * 0.05.  Other  parameters  are:**  7.3xl0~4,  0e/0j.  = 1, 
and  = 10-9.  The  numbers  on  the  contours  give  the  growth 
rate  in  units  of  2xl0-7wpe.  The  feature  on  the  left  is 
the  EID  instability  (hereFin  the  resonant  backscattering 
regime) . The  butterfly-like  pattern  centered  about  the 
point  kQ-k  = kof  is  the  SM  instability,  and  the  small 
structure  perpendicular  to  this  same  point  is  the  OTS  in- 
stability (here  in  the  regime  described  under  case  A in 
the  text).  These  parameters  are  appropriate  to  a pump 
generated  by  an  electron  stream  associated  with  a Type  III 
solar  radio  burst  at  1 (about  7xl09  km  from  the  sun). 
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90* 


o* 


Figure  2b.  Contours  of  constant  growth  rate,  Im  a , 
as  a function  of  the  wavevector  of  one  of  the  excited 
Langmuir  waves  (see  Figure  one).  Here  WQ  * 1.4x10-6  and 
k0  = 0.026.  The  other  parameters  are  the  same  as  in 
Figure  two  a) . Although  it  is  now  obscured  by  the  other 
contours,  about  the  point  - k * is  a structure 
qualitatively  similar  to  tKat  in  the  previous  figure, 
which  is  the  OTS  instability.  These  parameters  are  ap- 
propriate to  a pump  generated  by  an  electron  stream  of  the 
kind  associated  with  a Type  III  solar  burst  at  30  RQ. 


90* 


0* 


Figure  2c.  Contours  of  constant  growth  rate.  In  U, 
as  a function  of  the  wavevector  of  one  of  the  excited 
Langmuir  waves  (see  Figure  one).  Here  W = 5xl0~6,  and 
kQ  = 0.013.  The  other  parameters  are  th8  same  as  in 
Figure  two  a).  The  EID  instability  is  here  a resonant, 
forward  scattering  instability,  marked  by  contour  ”3" 
on  axis.  These  parameters  are  appropriate  to  a pump 
generated  by  an  electron  stream  of  the  kind  associated 
with  a Type  III  solar  radio  burst  at  215  RQ. 
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Figure  three.  Solutions  to  the  dispersion  relation, 
equation  (6),  showing  contours  of  constant  growth  rate, 

Im  «,  as  a function  of  the  wavevector  of  one  of  the  daughter 
Langmuir  waves,  kp—k.  These  solutions  all  occur  in  a 
regime  intermediate~between  cases  B and  D for  the  OTS  in- 
stability, where  10ko2  «.  W„  and  WQ  « 0(64w).  Here  WQ  = 0.04 
and  kQ  = 0.0067.  The  othe  ■ parameters  are  given  in  Table  V. 
These  parameters  are  those  of  a pump  generated  by  an  auroral 
streamer.  In  this  regime,  the  EID  has  disappeared  and  the 
the  OTS  instability  is  on  axis.  The  numbers  on  the  contours 
give  the  growth  rate  in  units  of  2.5x10-4  u>pe. 
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This  occurs  when  10ko2  is  on  the  order  of  WQ.  A set  of 
parameters  in  this  transitional  regime  generated  the  solu- 
tions shown  in  Figure  four.  In  this  regime,  the  condition 
that  y > x is  only  marginally  satisfied  on  axis.  Qualita- 
tively, the  transition  occurs  as  the  growth  rates  on  axis 
and  perpendicular  become  nearly  equal  and  growth  is  almost 
isotropic  in  k.  As  Figure  four  shows,  in  this  transitional 
regime,  nowhere  is  the  growth  its  maximum  value  of  Wc/8 
for  the  OTS  instability,  even  though  the  (now-decreased) 
maximum  does  occur  at  the  S given  by  (16a).  Equation  (16a) 
describes  an  ellipse  about  the  point  kQ-k  = k^. 

These  features  of  the  k » 0(kQ)  reqime  (with 
|wj  cj;  ) can  be  understood  using  the  approximate  expression 
for  6 at  maximum  growth  rate  (from  (16a)  to  first  order  in 


WQ): 


2.kl  Sin' t V' 


As  soon  as  kQ  = 0(4k),  the  sin^  term  degrades  the  growth 
rate  at  y = 90°.  The  increase  in  the  growth  rate  for  small 
sin^/.,  shown  in  (16)  , then  competes  with  the  decrease  in 
growth  rate  that  comes  from  the  violation  of  the  condition 
y ■>  x = 3kk0  (l-sin^Jt) 15 , that  occurs  when  sin2)6  is  small. 
This  regime  is  the  transition  among  all  four  regimes  con- 
sidered here. 

C) . |w| » u);t  k#»^.  This  limit  occurs  when  C4*/t<<Kq 
and  Wo«10ko2.  In  this  case,  we  have: 
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Figure  four.  Contours  of  constant  growth  rate,  Ira  <*>, 
as  a function  of  the  wavevector,  kj,  - k,  of  one  of  the 
daughter  Langmuir  waves  (see  Figure  one).  Here  WQ  °6.5xl0-4. 
6.5xl0“4  and  kQ,  shown  by  an  arrow,  has  magnitude,  0.028. 
Other  parameters  are:  9 « 7.3xl0“4,  0e/Q^  = l,Yi.»  10“9. 

The  numbers  on  the  contours  give  the  growth  rate  in  units 
of  3x10“5.  The  contour  labelled  "3"  is  the  maximum 
growth  rate  for  the  EID  instability,  which  here  is  in 
the  resonant  forwardscattering  regime.  The  four  smaller 
lobes  at  acos  (^*)  * 15°  are  the  SM  instability.  The 
OTS  instability  here  is  intermediate  between  regimes  A and 
C of  the  text,  and  its  maximum  growth  rate  is  shown  by  the 
"1.7"  contour. 
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D (k, u)  ~ ( U)  *-  Ua7  * i Va  <*j)/ojS 


so  that  (13)  becomes: 


M.y  =o. 


where  we  have  assumed  that  Im  y >>  v, . The  solution  to 

! *i  " 

• i l- 

(17),  we  have  found  nunerioaiiy , is  only  applicable  for 
very  large  powers,  WQ>  1.  However,  it  does  give  good 
qualitative  descriptions  for  Wo^.10“2.  The  solution  for 
CJ  is : 

s*[*  i - \ i + <xw<yny  s1  ] • 

Again,  we  differentiate  with  respect  to  and  k,  and  get, 
for  this  limit,  that  } * 90°  for  maximum  growth 
independent  of  k as  long  as  k « kQ.  The  value  of  k is 
determined  by  a cubic  equation  for  8 (or  k2)  which  comes 
from  the  differentiation  of  (17a)  at  f = 90°: 

«»  nr-wffl'i&m-iw)-*' 

where  S9  * -3k02/2.  The  asymptotic  growth  rate,  for 
5 »W  is  (from  (17a)): 


(17c) 


Figure  five  gives  a solution  to  equation  (6)  in  this  re- 
gime. For  these  parameters,  we  predict,  from  (17b)  and 
(17c),  that  k = 0.044  and  Im  lo  — 6xl0“2.  The  numerical 
results  give  k = 0.038  and,  for  the  growth  rate,  Im u - 2xl0-^, 


1 


8 3 


Figure  five.  Contours  of  constant  growth  rate,  Im  w, 


as  a function  of  the  wave  vector,  Jcq  - k,  of  one  of  the 
daughter  Langmuir  waves  (see  Figure  one).  Here  WQ  = 0.15 
and  k_  *=  0.05.  The  other  parameters  are;  d » 7.3x10*4, 

Yu  * 10“9  and  9e/f?i  ■ 1.  The  numbers  on  the  contours  aive 
the  growth  rates  in  units  of  10"3.  In  this  regime  the 
OTS  and  SM  instabilities  have  coalesced  and  the  EID  in- 
stability is  still  centered  near  the  condition  given  by 
perfect  frequency  matching. 
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(See  Table  IV  for  these  results.)  The  conditions  that  must 
obtain  on  kQ  and  WQ  for  self-consistency  are: 

W0  » 64*/r  , W 0 kj~  . 

D)  . ju>l  >>  , k y>  fto  • This  is  the  large  amplitude 

pump  region  for  the  OTS  instability,  when  we  require 
WQ  » 64oC,  wQ  »10ko2.  For  a large  WQ,  relative  to  kQ2 , 
as  is  shown  in  Figure  three  , equation  (17)  gives  the 
asymptotic  growth  rate,  but  now  we  have  growth  concentrated 

A A 

on  axis,  k*kQ  = 1.  The  solution  to  (17),  however,  has  no 
maxima,  but  only  a asymptotic  limit,  for  large  S , of: 


Notice  that  for  large  8 , that  is  k»  kQ,  that  the 
daughter  waves  set  up  a standing  Langmuir  wave  field.  For 
the  parameters  of  Figure  three,  for  example,  we  have  com- 
puted the  relative  amplitudes  and  phases  of  the  Stokes 
and  anti-Stokes  waves  for  the  OTS  instability,  and  they 
are  equal.  Since  the  frequency  of  each  wave  is  almost 
equal:  •*  u % w*  - <*>  , and  since  they  have  nearly 

equal  and  opposite  wavenumbers,  they  combine  to  form  a 
standing-wave.  This  has  important  effects,  one  of  which 
we  will  examine  in  an  application  of  the  fields  generated 
by  auroral  streams. 

The  lower  bound  on  Wc  for  this  regime  is  given  by 
the  condition  that  k >>  kQ  and  M»uA  which  give: 

» 

■ 

. I 

. ..  j 


TABLE  IV  (continued) 
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We  » 10  ho 
Wo  » -p  ■ 

The  transition  between  k » kQ  and  k«kQ  when  |**»l»  , is 

shown  in  Figure  six.  The  features  evident  in  Figure 
two  are  again  the  qualitative  determinants  for  this  case, 
similar  to  Figure  four  for  smaller  . 


ELECTRON-ION  DECAY  INSTABILITIES 

In  contrast  to  the  OTS  instabilities,  the  decay 
instabilities  depend  critically  on  the  resonance  of  the 
low  frequency  mode,  that  is:  Re  u>  * u/A.  Our  general 

schema  for  classifying  the  decay  instabilities  uses  this 
resonance  condition;  after  Fried,  et  al.  we  distinguish 
three  regions  of  k-space,  depending  on  their  resonance 
properties : 

1) .  resonant,  backscattering , when  kQ  > fz* /'j  , 

2) .  resonant,  forward  scattering,  when  ’fz»7}?kc>  f*7c  , 

3) .  ncnresonant,  forward  scattering,  k.0<Vw71T. 

Each  of  these  regions  must  be  examined  for  the  case  of 
the  low  frequency  resonance  being  a quasi-mode  (6e&  9^) 
and  for  the  case  of  a well-defined  mode  (9e  » 6^) . We  will 
examine  these  six  regimes  separately. 

Physically,  the  EID  instabilities  are 
the  decay  of  a plasma  wave  into  a plasma  wave  at  lower 
frequency  and  an  acoustic  mode  (or  quasi-mode) . However, 
to  treat  these  instabilities  analytically  is  not  entirely 
straightforward.  The  problem  is  one  of  approximating 
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Figure  six.  Solutions  to  the  dispersion  relation, 
equation  (6),  showing  contours  of  constant  orowth  rate,  Im  w , 
as  a function  of  the  wavevector  of  one  of  the  daughter 
Langmuir  waves  (see  Figure  one) . Here  W * 1x10"  and 
v , shown  as  an  arrow,  has  magnitude  0.032.  Other  para- 
meters are:  fc  = 7.3xl0"4, Vi.  * 10"9,  and  ©e/®i  = The 

numbers  on  the  contours  give  the  growth  rate  in  units  of 
2.5xi0"4.  In  this  large  WQ  regime,  the  OTS  and  SM  in- 
stabilities have  merged  ana  the  OTS  instability  itself 
is  almost  isotropic  in  k.  The  maxima  labelled  "3”  are 
the  EID  instability.  The  OTS  instability  is  in  a reqimc 
intermediate  among  the  four  cases  described  in  the  text. 
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D(k,  0)  for  a range  of  ion  temperatures.  When  ©j,  ^ , 

particle  effects  predominate  and  the  EID  instability  would 
be  more  properly  called  "induced  scattering  off  ions , " 2 , 11 , 15 
since  the  low  frequency  resonance  is  not  well  defined. 

When  0e  » 0^ , however,  the  acoustic  mode  is  well-defined. 

We  have  derived  an  approximation  which  treats  the  Ge^O; 
and  regimes  reasonably  well,  using  a recursive  approach 
due  to  DuBois  and  Goldman.*1  We  begin  by  noting  that  the 
LID  instabilities  obtain  when  only  one  high-frequency 
dielectric  function  is  near  resonance: 


+£a;)  >>  £ (ft -fci,  Uj  -U\) 

so  that  re  could  write  equation  (6)  as: 

(is)  eNL  (Is  k,'Lj-ue)^  — - c, 

if: 

k*6 « < , 

The  growth  rate  (or  damping)  of  this  new,  non-linear  mode 
is : 

y**-  ~ - 


k7  6 (k,w) 


The  maximum  qrov/th  rate  depends  on  the  maximum  of: 

[ k 6 ( k , cJ]  | 

We  will  treat  the  0e»  0^  case,  in  which  there  is  not  a well- 
defined  mode  at  the  acoustic  frequency  by  choosing  WA'  and 


I 


VA’  in  the  fluid  form: 


01 


(18a) 


in  such  a way  that  the  Im  k2€(k,  W)  “ ^ for  the  fluid  case 


aqrees  with  the  exact  result.  Now: 

iW  w)]  (iuid. 


■ max 


2V/ 


The  fluid  model  attempts  to  approximate  the  bell-shaped  k- 
depenuence  of  D(k,  uj  ) near  the  acoustic  frequency  with  a 
i road  Lorentzian.  The  exact  result  is: 


-I 

jtxac*  -vnax 


^ 0.53 


We  will  attempt  to  use  the  fluid  model  by  normalizing  the 
fluid  Lorentzian  so  that  the  maximum  value  of  (18a)  agrees 
with  the  exact  result,  and  so  that  this  maximum  occurs  at 
the  same  value  of  k;  this  requires: 

v/  = 0'90coa' 

1.70  • 


(18b) 


These  are  the  values  (in  units  of  «pe  and  k^)  which  we 
will  use  if  9 G*0i*  T^e  results  quoted  in  the  first 
chapter  can  be  used  when  With  this  approximation 

scheme  for  D(k,  u) , let  us  look  at  these  six  regimes  for 
the  electron-ion  decay  instabilities. 

A).  Resonant  backscattering , ©c>>9^.  for  this 
temperature  range,  the  acoustic  mode  is  well-defined, 

Re  &ja>>  X a = va/2,  and  we  can  use  a resonance  approximation 
for  D (k , £v)  : 


and  we  have : 


6 Wo  + ^)  *2  4 S , 

so  that  equation  (6)  becomes: 

(2D  (V'OOa* iyA)( W'Xo-? +»yL) 

1 6 

This  equation  has  been  studied  in  detail  by  DuBois  and 
Goldman,  5 and  we  note  here  only  the  simplest  of  their 
results.  The  maximum  growth  rate  occurs  when: 

00  - U)A  = X«,  ♦ S = 8_  ■ 

which  is  merely  the  frequency  matching  condition: 

We  (k»)  = W,  (k)  t U)L  (fak). 

For  a Langmuir  pump,  this  can  be  written: 

(22)  k = 2Ao  cos  / - * 

The  requirement  for  resonant  backscattering  is  then  a 
condition  on  kQ,  since  here  k*kQ  = +1: 

he  yft/T  ■ 

Note,  however,  that  in  addition  to  the  well-known,  solution 
with  cos  /•  =1,  that  (22)  actually  defines  a locus  in 

phase  space;  it  will  turn  out  that  only  the  cos  ^ = 1 
solution  is  truly  a decay  instability,  and  so  we  defer 
discussion  of  the  other  parts  of  the  locus  to  the  section 


on  the  SM  instabilities.  The  growth  rate  on  this  locus 


fL_  , . - _ (>a+Vl)  rwTyi*  ' U)A  , t 

(22a)  JnOO-  — • 

The  threshold  is,  then,  approximately: 

Wc>  “ 1 6 fa  / , 


and  the  maximum  growth  rate  is: 


/ K«*a,*«* , 


(^‘Z 


5I^W*(f) 


< y„«yL 


The  qualititive  features  of  the  resonant,  backscattering 
EID  instabiltty  are  evident  from  Figure  two  a) , although 
this  is  for  a different  temperature  ratio.  The  structure 
on  the  left  of  the  graph,  showing  growing  backscattered 
waves,  is  the  EID  instability.  Since  the  maximum  growth 
is  on  axis,  JU*  sryU+a*  and  the  detailed  angular  depend- 
ence of  yU  only  serves  to  determine  the  spread  of  the 
backscattered  region  of  growth. 

B) . Resonant  backscattering,  Te  ~ . In  this 

case,  the  resonant  approximation  for  D(k,  w)  is  not  ac- 
curate; to  get  approximate  analytic  results,  wc  use  a 
normalized  Lorentzian  form  for  D(k,  w)  given  by  (18a) 
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and  (18L).  Our  numerical  work  was  done,  however,  with  the 
exact  expression  for  D(k,  <o)  . Approximately,  then: 

l jL.1 


(24)  (w’-or/*  iu>v;Uu>-8  ■ 

& 


- o 


We  have  been  able  to  solve  this  equation  in  the  frequently 
applicable  limit  that  yL  is  negligible  (i.e.,  much  smaller 
than  either  Ira  to  or  yA ) . To  treat  (24),  we  look  at  two 
regimes,  Im  yA  , and  Im  to  » yA  . In  the  latter, 
high  power  regime,  (24)  reduces  to  equation  (21)  of  the 
previous  case.  In  the  low  power  limit,  a perturbation 
analysis,  in  the  small  parameter  y/j^  gives  the  follow- 
ing results  (again  we  restrict  ourselves  to  the  case  that 
/Vs  • which  is  easily  satisfied  a posteriori). 

To  first  order  in  y/yA  , the  imaginary  part  of  equation 

(24)  gives: 

(25)  Re  U>  * S_[l  - P )]  J 

ans  substitution  of  this  into  the  real  part  of  (24)  gives: 

n , . _ yu.u*1  r s-v± ] . 

(26)  w g [ (SKUiY  * K'S.2  1 

We  can  maximize  Im  Oi  with  respect  to  k,  and  we  derive  a 
condition  analogous  to  (22)  , a wavenumber  matching  condi- 
tion: 


(27) 


A * 2ko  COS  f “ ' 


From  this,  the  imaginary  part  of  can  be  calculated  and 


maximized,  when 


A’ 


1: 


95 


(28,  JU(o|m„=AW./7.5/  • 

Figure  two  a)  shows  these  results.  The  growth  rate  and  k 
for  maximum  orowth  are  predicted  very  well  by  (27)  and 
(28)  (within  5%) . The  condition  on  kQ  for  resonant  back- 
scatter  is  easily  gotten  from  (27)  : 

C)  . Resonant,  forward  scattering,  fle  >V  6^.  In 
the  backscattering  case  considered  above,  the  frequency 
matching  condition  resulted  in  the  condition  on  the  wave- 
numbers  that: 

k-2kcU)S  /■  - - ffiT  , 

and  kQ  - k < 0.  However,  we  can  still  have  resonant  scat- 
tering in  the  regime  where  kQ  - k >0,  or,  in  terms  of  kQ: 

(29  fn  > K ■ 

When  this  is  the  case,  the  daughter  Langmuir  wave  is  scat- 
tered in  the  forward  direction;  however,  the  instability  is 
still  described  by  (21)  and  the  growth  rates  given  by  (23). 
Notice  that  the  growth  rates  may  differ,  since  k is  smaller 
in  this  regime  and  the  growth  rates  depend  on  k.  Figure 
four  presents  numerical  solutions  to  the  full  dispersion 
relation,  equation  (6)  in  the  equal  temperature  regime. 

The  numerical  results  for  the  regime  here,  0e»0±r  are 
qualitatively  very  similar. 
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D).  Resonant,  forward  scattering,  0e  ~ 6^.  As 
for  Te » , the  analysis  from  the  resonant  backscattering 

for  equal  temperatures  can  be  applied  to  the  forward 
scattering  case  with  equal  temperatures.  The  condition 
that  kQ  must  satisfy  is: 


(29a) 


> kc  > 0451*} 3 • 


At  this  point,  we  must  examine  a feature  of  equation 


(7)  which  we  have  thus  far  ignored.  We  begin  by  noting 
an  important  property  of  the  symmetries  in  (6) ; if  we  make 
the  transformation: 

k — •-& 

Cu  ► - CO* 

the  equation  remains  unchanged.  That  is  for  a given  growth 
rate  at  (Re(w0-v),  , the  growth  rate  at  (Re  (wQ+w)  ,kQ+k) 

will  be  the  same.  In  Figure  four  this  symmetry  is  evident, 
and  the  EID  growth  is  present  for  wavenumbers  larger  and 
smaller  than  kQ.  Thus,  strictly  speaking,  in  Figure 
two  a),  for  example,  we  should  have  drawn  a forward  lobe 
corresponding  to  the  backscattered  instability,  gotten 
by  inversion  of  the  backscattering  lobe  throuqh  the  point 

ho"!i  = lSo* 

To  understand  why  this  other  decay  instability  only 
exists  under  certain  conditions  (whose  exact  nature  is 
closely  tied  up  with  the  SM  instability),  even  though  the 
symmetry  in  (C)  seems  to  demand  it  always,  we  must  go 
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Lack  to  equation  (1),  and  ask  the  question:  What  are  the 

amplitudes  of  the  excited  fields?  To  answer  this,  we  use 
a result  from  normal  mode  theory,  noting  that  the 
solution  of  (6),  our  dispersion  relation,  is  essentially 
a solution  for  the  eigenvalues  of  the  matrix  in  (1).  It 
is  well-known  that  this  normal  mode  analysis  can  predict 
the  relative  amplitudes  of  the  excited  waves,  as  well  as 
the  normal  mode  frequencies.  These  amplitudes  will  be 
given  by  the  norm  of  the  eigenvectors  associated  with  the 


zero  eigenvalue 

. That 

is, 

given 

a mati 

h‘ 

a u 

a» 

/ £l ' 

a* 

6r 

E, 

^a„ 

e>l 

\h 

- 0 


then  we  know  for  non-trivial  : 


__  ~ 6»  Six 

djld't-Cjdij 


For  the  matrix  in  (1),  we  have: 


£ (£-&>,  m-ojA  „ tCk+ko^rcoo) 

E ^ JU*  ~ £(&-&>  t to- UJ») 


As  either  a look  at  the  frequency  matching  condition  or 
the  numerical  results  in  Table  IV  show,  only  the  branch  of 
the  LID  instability  involving  the  Stokes  wave  is  near 
resonance,  and  for  that  case: 


(30a) 
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When  (30a)  is  satisfied,  only  the  LID  instability  growth 
rates  occurina  on  the  left-hand  side  of  our  graphs  will 
correspond  to  a wave  with  a large  parametrically  excited 
amplitude.  In  the  opposite  case,  when: 

|Rfl.  6 Ck'ko,(^-M,)  ~ Rt  €(fc+k*,UJ«-aA>)j<  ^mCCUi-CJo), 

then  both  the  Stokes  and  anti-Stokes  modes  may  simul- 
taneously be  resonant.  We  can  then  write  (using  (ID): 

(30b)  |S+  ♦ S_l  * W LO  ■ 

If  (30b)  is  satisfied,  then  our  three-wave  treatment 
breaks  down.  A treatment  of  the  case  when  (30b)  is  sati- 
fied  is  taken  up  in  the  next  section,  on  the  SM  instabi- 
lity. 

E).  Nonresonant,  forward  scattering.  As  kQ  de- 
creases further,  and  kn  <.  S fTec/T*  * it  becomes  impossible 
to  satisfy  the  resonance  condition  for  the  Langmuir  wave, 
but,  the  growth  rate  of  the  EID  instability  does  not  drop 
off  sharply  once  kQ  decreases  beyond  H . To  treat  the 
nonresonant  case  analytically  is  difficult,  since  the 
equation  we  must  then  solve  is  intrinsically  quartic. 
However,  the  growth  rates  (approximately) , and  a quali- 
tative  feeling  for  the  depence  of  the  instability  on  kQ, 
can  be  obtained  in  the  limit  that  the  dampings  of  the  waves 
are  much  less  than  the  growth  rates,  while  W0<*  « 1.  In 
this  somewhat  artificial  limit,  (6)  can  be  written: 
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(31)  (CO*-  US)  (OJ-X,)1  - - o 

Now,  if  Re  xQ,  then: 

(jWw)1  = iky£i£-(-  ) ■ 

As  ko  decreases,  Imu  decreases.  Figure  seven  shows  sol- 
utions to  equation  (6)  in  the  nonresonant  forward  scatter- 
ing case,  although  here,  k < kQ.  It  should  be  noted  that 
a recent  paper  by  Fried,  et  al.  treats  the  decay  insta- 
bilities driven  by  a Langmuir  pump  in  one  dimension. ^ They 
present  very  detailed  results  on  the  dependence  of  the 
growth  rate  on  kQ. 

STIMULATED  MODULATIONAL  INSTABILITIES 

Under  the  name  stimulated  modulational  instabili- 
ties, we  include  those  interactions  which  require  the 
Stokes  and  anti-Stokes  high  frequency  waves  and  a resonant 
low  frequency  response  be  treated.  In  this  case,  three 
dimensional  effects  are  critical*  in  fact,  the  instability 
cannot  exist  in  one  dimension.  (Notice  that  the  OTS  in- 
stability, while  it  involved  both  high  frequency  waves, 
had  a non-resonant  low  frequency  response.  The  OTS  is 
properly  called  a modulational  instability,  but  not  a 
stimulated  modulational  instability.)  Inclusion  of  a 
third  resonant  wave  complicates  an  analytic  treatment,  and 
we  have  been  able  to  treat  this  case  only  in  the  very  low 


Figure  seven.  Solutions  of  the  dispersion  relation, 
equation  (6),  showing  contours  of  constant  growth  rate,  Im  w. 
as  a function  of  the  wavevector,  - k,  of  one  of  the 
daughter  Langmuir  waves  (see  Figure  one).  Here  WQ  = 4xl0~6 
and  has  magnitude  0.005.  Other  parameters  are  the  same 
as  for  Figure  six.  The  numbers  on  the  contours  qive  the 
growth  rate  in  units  of  2xl0~7.  In  this  regime,  the  EID 
instability  is  no  longer  resonant,  and  hence  has  the  much 
degraded  maxima  at  the  two,  small  contours  labelled  "1." 

The  other  maxima  is  the  OTS  instability,  here  almost  iso- 
tropic in  k. 


and  very  high  power  limits. 
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Even  so,  the  general  result  of  the  inclusion  of  both 
the  Stokes  and  anti-Stokes  waves  is  a lowering  of  the 
growth  rate.  This  an  effect  well-known  in  nonlinear 
optics  and  was  noted  in  the  case  of  a plasma  instability 
by  DuBois  and  Goldman. 3 They  showed  that,  in  fact,  for 
an  electromagnetic  pump,  the  threshold  for  the  four-wave 
instability  is  higher  by  a factor  yA  . We  have  used 
perturbation  methods  to  study  the  effect  of  the  coupling 
of  both  high  frequency  waves  and  have  been  able  to  show, 
directly  that  in  the  limit  y<tyA  and  the  limit  y » that 
the  coupling  of  the  anti-Stokes  wave  reduces  the  growth 
rate  of  the  instability.  This  seemingly  general  result  is 
essential  in  understanding  the  SM  instabilities. 

We  have  gotten  analytic  results  for  these  four-wave, 
resonant  instabilities  in  two  regimes.  In  the  low  power 
regime,  when  we  have: 

& » ^ w >>  Yl 

Figures  two  a),  on  page  77,  and  four,  on  page  81 » give 
numerical  results.  In  the  high  power  limit,  when  u,"»  , 

the  SM  instability  reduces  to  the  high  power  limit  of  the 
OTS  instability  (case  C and  D) ; analytic  results  deduced 
for  that  case  also  describe  the  SM  instability  in  the  high 
power  limit.  Essentially,  as  Figure  five  (page  83)  shows, 
the  SM  and  OTS  instabilities  merge  for  large  powers. 
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A striking  feature  of  the  numerical  results  for  the 
lower  power  regime,  as  shown  in  Figures  two  a)  and  four,  is 
the  wing-like  structures  at  1^  - k E These  lobes  are 

the  SM  instability.  We  begin  with  equation  (13)  written 
as,  using  the  ordering  yL  < < cj  «.  yA  : 

(32)  *\CJV/k)(lAJ-S.)(oO*S,)*VJoSuAz/<1/^  -o 

(it  turns  out  that  even  off  axis  the  approximation  that 

1*1  = yU+2  is  adequate  to  second  order  in  <yfa  for  the  small 

W here).  The  solution  to  (32),  to  first  order  in  y,  is 
o 

given  by  equations  (25)  and  (26)  and  their  anti-Stokes 


counterparts : 


fefci-  ±sTfi  - . 

fv.  , . s.  _] 

^ ' <5  (V-WJ1*  J *' 


and  these  have  their  maximum  growth  at: 


where  we  have  defined: 


1 TJTj 


Notice  that  when  we  take  the  upper  sign  in  (34)  , we 
recover  the  three-wave  EID  results.  The  reason  for  this 
is  clear  if  we  notice  that,  to  first  order  in  y: 

(w- $-)(<*>♦§*)  *$<■ 

" 2 S [fcj  + ■ 
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which  is  putting  one  high  frequency  wave  off  resonance  and 
fixed  at  (kQ±k,  w iwD)  = 28  . That  is  to  say,  four- 
wave  effects  are  also  second  order  in  y (or  WQ) . 

We  were  able  before,  retaining  only  three-wave 
effects,  to  predict  the  maximum  growth  rate  along  the  locus 
given  by  (35),  and  this  is  also  the  maximum  growth  rate 
for  the  SM  instability.  However,  the  three-wave  treat- 
ment always  predicts  maximum  growth  for  k = 0.  Our  numeri- 
cal work  shows  that  the  maximum  growth  rate  occurs  for 
small,  but  non-zero,  k.  This  is  because,  as  k goes  to 
zero,  we  have,  from  (30b) : 

I St  ♦ S.|  “ 2S  — 0 , 

so  that  both  high  frequency  waves  can  be  near  resonance 
(and  so  we  cannot  ignore  the  anti-Stokes  wave) . As  the 
interactions  becomes  truly  four-wave,  the  threshold  goes 
up,  the  growth  rate  goes  down;  and,  we  get  the  maximum 
growth  rate  off-axis.  There  are  two  competing  effects 
here : 

1) .  /*  has  its  maximum  when  k **  0;  that  is,  when 

and  k^tk  are  parallel.  Thus,  as  k goes  to  zero,  the 

effect  is  to  increase  fA  and  increase  the  growth  rate; 

2)  But,  when  k gets  smaller,  the  resonances 
of  the  Stokes  and  anti-~tokes  waves  get  closer.  This 
means  that  as  k-*0,  the  anti-Stokes  modes  is  more  and 
more  strongly  coupled,  which  has  the  effect  of  lowering 
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the  growth  rate. 

It  is  the  competition  of  these  two  effects  which  determines 
the  k for  which  maximum  growth  occurs. 

To  calculate  this  k requires  a second  order  pertur- 
bation analysis  of  (32).  We  have  done  this,  beginning 
by  writinq  (32)  in  the  form: 

(37)  (w'-UVWUftjfe  *l«j)  (ZS  ♦I'u)  t W,Sfci,*/i/4  - 0, 

using  the  definition: 

(M  ~ t e ♦ t y • 

We  have  four  small  parameters: 

e,  y,  W*. 

since  the  last  is  necessary  to  treat  the  dependence  of 
/MZ  on  k: 

(38)  /Ax»l-fely»7l , f = C0S-'{\X)  . 

The  result  of  this  perturbation  analysis  is  that  the  maxi- 
mum growth  rate,  on  the  locus  given  by  (35)  occurs  at  the 
k for  which: 

2S--S,  • 

(39) 

This  result  is  independent  of  WQ  for  low  powers.  It  pre- 
dicts the  k values  for  Figures  two  a)  and  four  within  5%. 

In  general,  the  locus  given  by  (35)  defines  an 
angle  at  which  the  SM  instability  will  occur  for  very  small 
k,  like  those  prescribed  by  (39),  namely: 

(40)  tos  f = 2^7  [: 

i — 
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This  is  also  a condition  on  k : 

o 


ko  > 


0.46  yZ&Th  , 0e  « 0,  , 


0.50  VzaTT  , 0e  » 0.  * 

It  was  these  conditions  that  assured  (see  (29)  or  (29a)) 
in  the  case  of  the  LID  instability,  that  the  instability 
was  a resonant  one.  For  kQ  less  than  this  critical  kQ, 
the  SM  instability,  now  on  axis,  becomes  part  of  the 
forward  scattering  nonresonant  EID  instability.  Thus,  for 
kc<  0 45\fz*V3  the  SM  and  EID  become  identical  instabilities, 
while  for  W0»  64P/T,  the  SM  and  OTS  merge. 

Figures  eiqht  through  ten  summarize  the  above 
results,  by  labelling  the  dominate  instability  for  a 
Langmuir  spectrum  as  a function  of  kQ,  Wq,  and  the  ambient 
plasma  temperature  ratio.  These  diagrams  give  qualitative 
results,  based  on  the  instability  with  the  largest  growth 
rate.  Two  results  are  striking: 

1)  . For  small  enough  kQ,  only  the  OTS  instability 
is  strong; 

2) .  For  large  ©e/f?i,  the  EID  instability  dominates 

all  other  parametric  processes.  These  results  are  only 
applicable  to  the  case  when  we  can  ignore  the  damping  of 
the  high  frequency  modes.  For  large  different  results 

are  possible  in  a few  regimes. 

APPLICATIONS  TO  BEAM-GENERATED  SPECTRA 


We  now  turn  to  several  applications  of  the 
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Figure  eight.  Qualitative  features  of  Langmuir 
pump  decay,  as  a function  of  WQ,  the  normalized  pump 
power,  and  k , the  pump  wavenumber.  Here  9„/e i = 1 and 
we  have  a hydrogen  plasma.  There  is  assumea  to  be 
negligible  damping  of  the  high  frequency  waves.  k0-k 
is  the  wavenumber  o f the  fastest  growing  daughter  wave. 
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Figure  nine.  Qualitative  features  of  Langmuir  pump 
decay,  as  a function  of  kQ  and  9 e/0 i.  Here  WQ  « 10“®  and 
we  have  a hydrogen  plasma,  which  we  assume  to  have  neglig- 
ible damping  of  the  high  frequency  waves. 
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dispersion  relation,  equation  (6),  which  describes  the 
linear  stability  of  a Langmuir  pump.  We  will  examine  three 
physical  situations  in  which  an  enhanced  spectrum  of 
longitudinally  polarized  waves  is  created  by  an  electron 
beam : 

1) .  Type  III  solar  radio  bursts.  These  are  radio 
signals  associated  with  electron  streams  ejected  from  the 
sun  which  generate  electromagnetic  radiation  through  an 
intermediate  step  involving  an  enhanced  Langmuir  spectrum. 

2) .  Auroral  arcs.  The  ionospheric  electron 
streams  associated  with  these  arcs  create  a non-equilibrium 
Langmuir  spectrum  which  probably  acts  as  a pump  for 
parametric  processes. 

3) .  Laboratorv  electron  beams.  In 
experiments  using  a non-relativistic  electron  beam 
travelling  through  an  almost  homogeneous  plasma  with  no 
magnetic  field,  parametric  effects  due  to  beam  generated 
Langmuir  spectra  are  sometimes  evident. 

One  of  the  commonest  sources  of  enhanced  Langmuir 
radiation  is  an  electon  beam  travelling  through  a plasma. 
Through  quasilinear  growth  and  stabilization, this 
beam  generates  asspectrum  with  characteristics : 

(41)  w , Ako/k,  « AUw/Vi. , 

(42)  Vs/  < Jli-  fo- 

* * Vi* 

where  vb  is  the  beam  velocity,  Avb  the  spread  in  the  beam, 


110 


AkQ  is  the  spread  in  the  spectrum,  and  nj-,/ne  is  the 
ratio  of  beam  density  to  ambient  plasma  density.  Equation 
(41)  is  a statement  of  the  fact  that  the  beam  generates 
waves  through  a "resonant"  wave-particle  interaction, 
which  requires  the  phase  velocity  of  the  waves  to  be  the 
same  as  the  velocity  of  the  particles.  Equation  (42) 
is  an  upper  bound  on  the  energy  that  the  beam  can  trans- 
fer to  the  wave  spectrum. 

Strictly  speaking,  the  dispersion  relation  of 
these  beam  generated  waves  is  not  Langmuir-like;  accord- 
ing to  O'Neil  and  Malmberg , the  dispersion  relation  of 
these  "beam  modes"  is  given  by  a product  of  the  plasma 
dispersion  relation  and  a factor  due  to  the  beam:^0 

-KVk*  ikAVk)1  - Hk/ne  • 


They  have  studied  this  equation  in  detail,  and  their  re- 
sults for  weak  beams  may  be  summarized  by  noting  that  when 
the  following  inequality  is  satisfied: 

I St  < c 

(z *./  V.  V, 

then  we  can  conclude  that  the  dispersion  relation  of  the 
waves  is  essentially  that  of  Langmuir  waves,  with: 

u «■  (;  +3kS)'/l  «■  (?("*/«.)• 


This  condition  is  satisfied  for  all  of  our  applications 
except  the  last.  Since  it  is  also  the  case  that: 


) 


Ill 


we  have  ignored  beam  corrections  to  the  Langmuir  dispersion 
relation  of  the  pump  in  this  limit. 

We  are  making  a further  approximation  in  applying 
a single-mode  pump  theory,  when  we  apply  our  dispersion 
relation,  equation  (6),  to  the  case  of  a finite  bandwidth 
pump.  Mathematically,  we  are  making  an  assumption  about 
the  narrowness  of  the  Langmuir  spectrum  in  k-space;  v/e 
are  requiring  essentially,  that  the  spectrum  be  narrower 
in  k than  the  high  frequency  dielectric  function  (see 
equation  (9)  for  an  exact  version  of  this).  The  effect 
of  a broad  pump  in  three  dimensions  is  discussed  in 
reference  8.  Here,  we  only  note  the  most  important  result, 
namely  that  the  condition  on  the  spectrum  can  be  written: 

|w  T Ik7  -3k  k,l  >>  J 3k  Mj  , 

where  it  is  important  to  realize  that  it  is  the  spread  of 
the  spectrum  in  the  direction  of  the  fastest  growing  mode 
that  appears  in  this  equation.  Note  also  that  the  left 
side  of  this  inequality  depends  on  the  growth  rate  of  the 
instability . 

In  the  further  approximation  that  the  beam-wave  in- 
teraction is  unimportant,  we  can  study  the  properties  of 
the  spectrum  it-self.  This  will  be  appropriate  in  a number 
of  circumstances.  If  for  example,  the  beam  passes  through 
the  system  quickly,  but  the  Langmuir  spectrum  remains 
at  a high  level  for  times  long  compared  to  the  beam  duration 
or,  if  the  beam-wave  interaction  has  characteristic  times 
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much  lonqcr  than  the  growth  rates  predicted  by  equation  (C), 
then  the  predominant  processes  in  the  beam  evolution  will 
be  parametric.  In  any  case,  given  an  assumed,  enhanced 
spectrum,  equation  (6)  gives  growth  rates  and  wavenumber 
shifts  for  the  linear  relaxation  of  the  spectrum. 

While  our  dispersion  relation,  equation  (6) , will 
describe  the  simplest,  linear  relaxation  of  the  beam-en- 
hanced spectrum,  our  formalism  does  not  directly  give  any 
insight  into  the  actual  evolution  of  the  beam-plasma  system. 
Our  methodology  in  the  following  applications  is  to  use 
the  properties  of  the  daughter  waves  predicted  by  equation 
(6)  to  make  qualitative  statements  about  the  importance 
of  linear,  parametric  effects  on  the  beam-plasma  evolution. 
We  can  examine,  for  example,  whether  spectral  evolution 
might  be  able  to  compete  with  quasilinear  relaxation  of  the 
beam,  if  it  turns  out  that  the  parametric  qrowth  rate  is 
comparable  to  the  quasilinear  growth  rate.  Our  theory  of 
the  linear  relaxation  of  the  spectrum  offers  important 
qualitative  insights  into  features  of  these  beam-plasma 
interactions,  but  it  cannot  make  quantitative  predictions 
of  the  detailed  evolution  of  the  beam  or  spectrum.  Thus, 
our  descriptions  of  these  three  applications  only  touches 
on  the  highly  nonlinear  processes  that  actually  determine 
the  physics;  but,  we  apply  (6)  in  an  effort  to  find  regimes 
where  parametric  effects  most  likely  play  a role  in  the 
total  (nonlinear)  interaction. 


We  take  three  cases  of  current  interest  and  compare 
the  results  of  numerical  and  analytic  solutions  of  (6)  with 
experimental  data  from  the  beam-generated  spectra  associated 
with  Type  III  solar  radio  bursts,  with  auroral  arcs,  and 
with  nonrelativistic  electron  beams  (with  no  magnetic 
field) . Table  V summarizes  the  physical  parameters  of 
each  of  these  beam  generated  spectra. 

Type  III  Solar  Radio  Bursts8'14'15'21-25 

It  is  believed  that  the  radio  discharges  of  solar 
origin,  at  200  mH  to  lOmH,  with  durations  of  one  to  ten 
seconds  and  having  a distinctive  harmonic  structure — Type 
III  solar  radio  bursts — are  due  to  high  energy  electron 
bursts  ejected  from  the  sun.  These  streams  of  electrons 
are  then  thought  to  interact  quasilinearly  with  the  sur- 
rounding plasma  (either  the  solar  corona,  solar  wind,  or 
interplanetary  plasma)  to  generate  an  enhanced  Langmuir 
spectrum,  which,  in  turn,  is  responsible  for 
the  observed  radio  signals.  The  role  of  the  parametric 
instabilities  in  the  evolution  of  this  beam-spectrum 
system  has  been  considered  by  several  authors,  the  EID 
instability  by  Zheleznyakov  and  Zaitsev15  and  the  OTS  in- 
stability by  Papadopoulos , et  al.22  and  by  Smith,  et  al.26 
The  basic  attack  on  the  problem,  from  a suggestion  by 
Kaplan  and  Tsytovich,14  is  to  study  the  relaxation  of 
the  Langmuir  spectrum  and  determine  whether  this  relaxation 
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occurs  in  such  a way  as  to  affect  the  quasilinear  decay  of 
the  beam.  Two  considerations  are  relevant: 

1)  . Do  the  parametric  processes  involved  in  the 
relaxation  of  the  Langmuir  spectrum  take  the  energy  of 

the  spectrum  out  of  resonance  with  the  beam  (and  so  destroy 
the  quasilinear  beam-spectrum  interaction);  that  is,  are 
either  of  the  following  inequalities  satisfied: 

Vl  - AVk  > V^‘Vk  ; 

(43)  -yk  ? 

2) .  Does  this  transfer  energy  occur  rapidly  enough 
to  compete  with  the  quasilinear  processes? 

The  second  of  these  questions  has  been  answered  partially 
by  references  8 and  15,  where  it  is  shown  that  neither  the 
OTS  or  EID  instability  occur  fast  enough  to  actually 
prevent  quasilinear  decay  of  the  beam. 

With  our  dispersion  relation,  (6) , we  have  studied  the 
behavior  of  the  enhanced  Langmuir  spectrum  and  have  found 
some  parameter  regimes  in  which  the  parametric  processes  can 
play  an  important  role  in  the  beam  dynamics,  even  though 
they  may  not  be  able  to  actually  prevent  decay  of  the  beam. 

We  begin  with  a numerical  solution  to  the  beam-spectrum 
interaction  in  the  case  of  solar  bursts,  recently  done  by 
Ma9blssen.27  This  numerical  work  solves  the  quasilinear 
equations  in  two  dimensions,  but  assumes  a one-dimensional 
spectrum  and  neglects  the  nonlinear  relaxation  of  the 
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spectrum.  From  Magelssen's  data  we  have  taken  the  para- 
meters necessary  for  a study  of  the  spectrum,  namely,  kQ 
and  WQ,  as  a function  of  the  distance  from  the  sun  that  the 
electron  pulse  has  travelled  (after  it  has  left  the 
surface  of  the  sun) . Table  V includes  the  parameters  for 

rn 

three  distances  from  the  sun,  and  Figure  two  (pages  57  and 
58)  show? contour  plots  of  the  parametric  growth  rates 
predicted  by  equation  (6).  We  have  used  Magelson's  data 
rather  than  the  approximations  from  a homogeneous  theory 
of  the  beam  generated  waves  given  by  (41)  and  (42)  for 
kQ  and  WQ. 

Magelssen's  results  show  that  kc  is  about  0.05 
(in  units  of  the  Debye  wavenumber)  near  the  sun,  and  kQ 
decreases  to  0.016  at  the  earth's  orbit  (about  215Re, 
where  1RC  is  one  solar  radius,  a distance  of  about  7x10^  km). 
His  calculations  of  the  Langmuir  wave  spectrum  intensity, 
which  include  effects  due  to  wave  reabsorption  and  "pile- 
up  give  typical  values  of  WQ  of  10  ° to  10  (the 
wave  intensity  normalized  to  the  thermal  enerqy  of  the 
plasma) . 

The  OTS  instability,  in  this  regime,  always  has  its 
maximum  growth  for  a k perpendicular  to  ^ and  for  values 
of  k much  less  than  kc.  Indeed,  it  would  require  stream 
densities  higher  by  several  orders  of  magnitude  to  result 
in  the  OTS  instability  having  any  growth  on  axis  for  even 
the  smallest  kQ  given  by  Magelsscn.  This  contradicts  the 
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results  of  both  references  22  and  26,  both  of  which  use 
the  OTS  instability  to  provide  for  energy  transfers 
which  satisfy  (4 3) and  so  could  account  for  beam  stabiliz- 
ation. We  believe  that  the  OTS  instability  will  have  a 
minimal  effect  on  the  beam-spectrum  system,  because  of  the 
impossibility  of  satisfying  (43)  for  the  powers  of  the 
spectrum  appropriate  to  Type  III  solar  bursts. 

However,  as  Zheleznyakov  and  Zaitsev  noted,  the 
EID  instability  can  transfer  energy  to  Langmuir  waves  with 
velocities  sufficiently  different  from  the  velocity  of  the 
beam  to  satisfy  (43)  only  when  it  is  not  a backscattering 
instability.  For  the  values  Magelssen  gives,  we  find  that, 
the  EID  instability  from  1 Rg  to  80  is  a backscattering 
instability.  From  80  RQ  to  450  Re,  somewhat  beyond  the 
earth,  the  EID  instability  is  a forward  scattering  instabil- 
ity. And  beyond  450  R , k is  so  small  that  the  condition 

o o 

obtained  from  combining  (43)  and  (28a)  is  violated  and  the 
EID  instability  becomes  a nonresonant  forward  scattering 
process.  Note  that  for  the  low  values  of  Wq  for  these 
spectra,  that  there  is  no  possibility  of  satisfying  the 
four-wave  coupling  condition,  equation  (20).  Thus,  if  we 
accept  Zaitsev  and  Zhelesnyakov' s argument  that  multiple 
scattering,  from  the  backscattering  EID,  prevents  its 
being  able  to  transfer  energy  to  nonresonant  waves  (that 
is,  waves  not  in  "resonance"  with  the  beam),  we  conclude 
that  only  between  80  and  450  can  the  EID  instability 
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shift  energy  out  of  resonance  with  the  beam.  We  also  note 
that  the  SM  instability  occurs  for  all  these  sets  of  para- 
meters for  distances  less  than  about  450  R , but  at  such 

© 

small  k-shifts  as  to  violate  (43). 

Our  conclusions  here  are  dependent  on  two  striking 
features  of  Magelssen's  results: 

1) .  The  small  values  of  WQ.  Many  authors  have 
conjectured  much  larger  values  of  WQ.  If  this  were  the 
case,  the  OTS  instability,  for  example,  could  have 
growth  on  axis,  if:  WQ  » 10ko2.  In  addition,  at  higher 
powers,  the  OTS  instability  will  result  in  an  angular 
spreading  of  the  spectrum,  even  if  kQ  is  not  small  enough 
to  satisfy  the  above  inequality  (this  effect  is  well 
shown  in  Figure  five,  page  83).  However,  to  get  into 
either  of  these  regimes,  even  for  the  smallest  of  the 

kQ  appropriate  to  the  solar  problem,  kc  = 0.013  (in  units 
of  kDe) , requires  a power  several  orders  of  magnitude  great- 
er than  Magelssen  finds.  If  WQ  were  on  the  order  of 
0.1,  it  is  even  possible  to  use  the  dipole  approximation 
for  this  problem,  as  Papadopoulos  has  done.22 

2) .  The  large  values  of  kQ.  Most  recently.  Smith, 
et  al.26  have  presented  a calculation  using  a value  of  k0 
of  0.013.  In  this  regime,  the  SM  instability  has  merged 
with  the  now  four-wave  EID  instability.  While  it  seems 
that  Smith,  et  al.  have  called  this  nonresonant  forward 
scattering  an  OTS  instability,  their  numerical  results 
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show  the  interaction  to  important  in  the  evolution  of 
the  beam  and  spectrum.  They  show  numerically  that  this 
growth  rate  shifts  energy  out  of  resonance  with  the  beam 
quickly  enough  to  stabilize  the  electron  beam.  Maqelssen's 
data,  however,  shows  that  the  value  of  k_  that  thev  use, 
as  appropriate  to  the  spectra  that  should  exist  near 
the  earth,  does  not  in  actuality  occur  except  for  distances 
on  the  order  of  twice  the  earth's  orbit  (near  450  R ). 

We  have  used  Magelssen's  data,  in  spite  of  these 
other  considerations,  because  his  results  so  closely 
mirror  the  best  observations  at  this  point;  the  essential 
effects  that  he  includes  are  the  inhomogeneity  of  the 
ambient  plasma,  and  the  finite  size  of  the  beam.  The  com- 
bination of  a sharp  beam  front  and  an  inhomogeneous  plasma 
can  stabilize  the  beam  without  the  inclusion  of  mode- 
couplinq  effects.  Even  with  his  very  low  values  of  WQ, 
the  observed  shape  of  the  spectra  are  nicely  accounted  for. 
What  is  necessary  now  to  improve  his  numerical  work,  is 
the  inclusion  of  these  mode-coupling  effects.  Our  results 
above  give  some  guide  to  the  regimes  where  these  effects 
will  be  most  important. 

While  our  qualitative  considerations  seem  to  show 
that  the  parametric  instabilities  cannot  stabilize  the 
electron  stream  associated  with  Type  III  bursts  (at  least 
between  the  earth  and  the  sun) , it  has  been  theorized  that 
they  Le  important  in  understanding  the  production  of 
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electromagnetic  radiation.  The  critical  consideration 
ie  the  conversion  of  Langmuir  radiation  to  electromagnetic 
radiation,  at  the  second  harmonic,  for  example  (the  char- 
acteristic feature  of  Type  III  radiation  is  the  presence 
of  both  the  fundamental  and  second  harmonic) , is  the 
angular  relations  between  the  Langmuir  wavenumbers  which 
can  allow  coalescence  of  two  Langmuir  waves  to  form  a 
transverse,  electromagnetic  wave  at  twice  the  plasma  fre- 
quency. (The  electromagnetic  ratiation  at  the  plasma 
frequency  is  usually  attributed  to  scattering  by  ion  flue- 
tuations.  J We  will  not  consider  that  process  here.) 
Conservation  of  energy  and  momentum  in  this  wave-wave 
process  require  that  the  wavenumbers  of  the  Langmuir  waves 
be  approximately  equal  and  opposite.  Our  results  for  the 
OTS  and  SM  instabilities  in  this  regime  show  that  this  is 
impossible  to  obtain  for  the  daughter  waves  (contrary  to 
the  findings  of  Papadopoulos,  et  al.22).  In  addition,  our 
results  show  that  the  EID  instability  is  only  a back- 
scattering  process  in  which  kQ-k  * -kQ  for  distances  with- 
in 16  R0.  Beyond  this  point,  kQ < 0.03,  and  we  have,  from 
equation  (27) : 

k„-  k • - k,  *■  0. 9 il*/ 3 Z-O.OI  . 

making  kQ-k  and  kQ  unsuitable  for  this  process  of  generat- 
ing radiation  at  26^e.  Our  analysis  shows  that  the  observed 
radiation  at  the  second  harmonic  must  be  due  to  processes 
other  than  parametric  ones. 


Auroral  Arcs}  28, 29 

Papadopoulos  and  Coffey  have  suggested  the  OTS  in- 
stability as  a mechanism  involved  in  the  production  of 
another  astrogeophysical  phenomenon,  the  auroral  arcs,  in 
a situation  similar  to  the  Type  III  solar  bursts,  experi- 
mental data  show  distinctive  electron  distributions  pro- 
pagating along  magnetic  field  lines  in  the  ionosphere, 
which  are  associated  with  auroral  discharges.  There  are 
two  pieces  of  data  that  naive,  homogeneous,  quasilinear 
theory  does  not  predict:  the  relative  stability  of  these 

beams  and  the  anomalous  resistivity  that  they  encounter 
along  the  magnetic  field  lines. 

Our  analysis  of  the  linear  stability  of  a spectrum 
of  Langmuir  waves  enhanced  by  a stream  typical  of  those 
associated  with  auroral  discharges  shows  that,  indeed, 
the  OTS  instability  will  have  a profound  effect  on  the 
resonance  relations  of  the  electron  beam  and  the  Langmuir 
spectrum.  Figure  three  (page  77)  shows  the  regions 
of  highest  instability  for  parameters  associated  with  an 
auroral  discharge  and,  given  the  fact  that,  k » kQ,  we 
have  a geometry  with  the  fastest  growing  mode  parallel  to 
kjo  . The  large  value  of  k allows  us  to  satisfy  (43),  even 
though  the  beam,  and  hence  spectrum,  are  relatively  broad; 
it  appears  that  the  OTS  instability,  in  one  dimension,  for 
these  parameters,  can  transfer  energy  from  regions  of  the 
Langmuir  spectrum  in  resonance  with  the  beam  to  regions  not 


in  resonance,  and  thss  possibly  stabilize  the  beam. 

However,  this  is  only  a necessary  condition  for 
stabilization.  As  Kaplan  and  Tsytovich  were  the  first  to 
point  out,  these  nonlinear  processes  also  must  proceed 
quickly  enouqh  to  change  the  spectrum  before  the  beam  has 
relaxed.  Zheleznyakov  and  Zaitsev  give  the  following 
condition  on  the  ratio  of  the  rates  of  quasilinear  relax- 
ation and  wave  transfer:^ 


(44)  e /r\  [c  We{i)  /w#a*o)  ] « I , 6 - *a  L '/V,.r. 

Using  the  results  for  the  wave  transfer  rate  from  equation 
(17b)  and  the  usual  quasilinear  relaxation  rat^,  we  get 
the  condition: 


(45) 


There  is  a discrepancy  between  (45)  and  the  results  of 
Papadopoulos  and  Coffey,  which  appears  to  come  from  their 
expression  for  the  maximum  growth  rate,  taken  from  Nishi- 
kawa,  which  is  incorrect  in  this  regime.  Although  data  on 
the  auroral  plasma  is  approximate,  we  know  that  A v^  /vb  is 
on  the  order  of  1/3  and  that  nb  * 0.02  cm"3,  so  that, 

(45)  gives:  ne  > 2xl05  cm~^.  This  is  100  times  greater 

than  the  result  of  Papadopoulos  and  Coffey,  and  an  exami- 
nation of  ionospheric  electron  densities  shows  that  this 
condition  is  never  satisfied  in  the  ionosphere.  The 
highest  electron  densities  observed  3®  are  about  10^,  which 


occur  at  about  300km,  much  higher  than  the  observed  height 
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of  auroral  discharges.  Thus,  we  conclude  that  are  no 
parametric  processes  sufficiently  fast  to  stabilize  the 
beams  associated  with  auroral  discharges.  However,  the 
other  important  features  of  these  discharges,  the  high 
anomalous  resistivity  and  power  law  tails  on  the  electron 
distributions,  for  example,  may  require  the  details  of  the 
linear  stabiltiy  of  a Langmuir  spectrum  to  be  fully  under- 
stood (and  Papadopoulos  and  Coffey  have  examined  a number 
of  these  questions),  but  in  a situation  where  the  OTS  in- 
stability does  not  stabilize  the  beam.  In  recent  paper, 
Matthews,  Pongratz,  and  Papadopoulos  correct  their  formula 
for  the  growth  rate  and  derive  a formula  similar  to  (45) 
for  nfa/ne.  They  also  report  new  data  which  shows  that 
ne  » 2xl05. 

We  have  developed  a formalism  sufficient  to  treat 
another  aspect  of  the  problem.  The  dynamics  of  the  Lang- 
muir spectrum,  after  an  initial  OTS  instability  transfers 
energy  to  large  wavenumbers,  kD  ± k x t0.05,  will  be 
determined  by  the  instabilities  of  these  daughter  waves. 

The  interesting  feature  of  this  situation  is  that  the  large 
wavenumber  daughter  waves  form  a standing  wave  field.  It 
is  straightforward  to  write  down  the  5 x 5 determinental 
condition  that  describes  the  coupling  of  the  low-frequency 
response  at  (k,  u)  and  the  high  frequency  waves  at 
(kQ  t k,  CO,t  U)  ) and  (kQ  ± k,  U/eTco).  The  kQ  is  now  the 
fastest  growing  mode  in  Figure  three,  and  the  WQ  will 
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be  determined  by  the  saturation  of  the  OTS  instability  in 
the  auroral  streamer.  On  this  value  of  WQ,  we  can  only 
speculate;  the  studies  of  the  nonlinear  saturation  of  the 
OTS  instability,  even  in  an  homogeneous  plasma  (which  the 
auroral  plasma  is  not)  are  vague. 30, 31  we  have  looked  at 
a range  of  values  of  WQ  from  1.6x10-6  to  4xl0~2  (figvred 
on  the  local  amplitude  of  the  field,  that  is,  the  sum  of 
the  amplitudes  of  both  components  of  the  standing  wave) . 

At  low  powers,  our  numerical  results  are  straight- 
forward in  that  we  get  directly  the  travelling-wave 
results  for  a pump  of  one-quarter  the  power.  For  the  kQ 
appropriate  to  the  daughter  waves  of  the  auroral  streamer, 

k =0.05,  the  numerical  work  we  have  done  for  a travel- 
o 

ling-wave  pump  of  the  same  kQ  (see  Figure  two  a),  pane  77) 
gives  the  same  results  that  would  obtain  for  a standing- 
wave  pump  with  WQ  = 1.6x10“^.  That  is,  as  we  saw  was 
possible  in  Chapter  II,  only  one  component  of  the  standing- 
wave  pump  is  driving  the  instability.  This  means,  most  im- 
portantly, that  the  auroral  streamer  parameters  generate 
a set  of  daughter  waves  which  must  be  treated  in  three 
dimensions,  even  though  the  parent  instability  occurs  all 
in  one  dimension.  The  relaxation  of  the  spectrum  of  the 
auroral  stream  occurs  in  three  dimensions. 

At  higher  powers,  WQ  > lxl 0“ 3 , we  have  been  unable 
to  solve  the  dispersion  relation.  The  exact  reason  for 
this  is  unclear,  but  there  is  evidence  that  our  expansion 
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in  I.Q  breaks  down.  The  major  difficulty  is  that  the 
celation  of  third  and  fourth  order  terms  in  CQ,  which 
occured  in  derivinq  our  travelling-wave  dispersion  relation, 
equation  (6),  does  not  occur  for  the  standing-wave  case. 

This  results  in  terms  remaining  in  the  dispersion  relation 
which,  as  WQ  gets  larqer,  begin  to  dominate  some  second 
order  terms.  The  conclusion  we  drew  above  about  the  im- 
portance of  three  dimensional  effects  does  not  seem  to 
be  affected  by  this  difficulty  in  our  formalism,  but,  for 
WQ  greater  than  lxlO-3,  a more  careful  study  of  the  para- 
metric processes  driven  by  a standing -wave  pump  must  be 
done . 

v>  34 

Laboratory  Electron  Beams  ' 

Quon,  et  al.  have  done  a series  of  laboratory  ex- 
periments involving  nonrelativistic  electron  beams  in 
which  parametric  processes  in  the  beam-produced  wave 
spectrum  are  evident.  Their  overall  results  show  that  the 
specific  parametric  interaction  depends  sensitively  on  the 
velocity  of  the  electron  beam,  or  equivalently,  on  the 
characteristic  wavenumber  of  the  enhanced  spectrum  of  plasma 
waves  that  the  beam  drives  up.  The  linear  stages  of  the 
decay  of  the  Langmuir  spectrum  that  this  beam  produces 
should  be  described  by  the  results  we  have  derived. 

We  first  note  that  the  parameters  of  Quon's  plasma 
are  very  different  from  the  astrophysical  plasmas;  the  mass 
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ratio  is  much  smaller  for  the  laboratory  situation  anc 
the  electron  temperature  is  much  qreater  than  the  ion 
temperature.  The  general  results  we  have  derived  for  this 
regime  show  that  the  LID  backscattering  instability  will 
dominate,  when  parametric  effects  arc  evident. 

In  the  experiments  done  by  Quon,  et  al.  two  bean 
velocities  are  examined.  In  the  case  of  the  slower  bean 
our  analysis  can  be  used  directly;  for  the  faster  beam, 
hwwever,  some  modifications  are  required.  In  the  case  of 
the  first  beam,  the  growing  solutions  to  equation  (6)  are 
shown  in  Figure  eleven  a) . The  large  wavenumbers  in  this 
regime  required  our  using  Landau  damping  for  the  high 
frequency  waves.  As  Quon  suggests,  the  dominate  process 
is  the  well-known  decay  instability,  and  the  parametric 
part  of  the  results  he  presents  are  described  by  the 
well-known  theory.  We  note  that  it  is  only  under  conditions 
of  small  ot  and  T»  I that  the  usual  results  (of  a back- 
scattered  wave  at  -ICq  and  ion  fluctuations  at  k^  = 2k0) 
obtain.  The  maximum  growth  rate  that  we  predict  from  our 
analytic  results  (equation  (23))  is  1.6xl0-3,  which  cor- 
responds to  the  lowest  value  of  WQ  that  Qron,  et  al . look 
at,  with  a beam  density  of  5%  of  background.  We  have  in- 
troduced no  phenomenological  constants:  the  damping  is 
purely  Landau  damping;  the  fields  have  been  calculated  from 
(41).  In  addition  to  the  decay  region,  our  numerical 
results  predict  another  region  of  growth.  Thi3  corresponds 
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Figure  eleven  a) . Solutions  to  the  dispersion 
relation,  equation  (6),  showing  contours  of  constant 
growth  rate,  Im  U,  as  a function  of  the  wavenumber  of 
one  of  the  daughter  Langmuir  waves,  kQ-k.  Here,  W = 0.0C; 
kQ  = 0.25.  Other  parameters  are;  o(  = 9x10"®,  0e/0?  = 15, 
and  the  high  frequency  deunping  is  the  sum  of  Landau  damp- 
ing and  a collisional  contribution  equal  to  3xl0~6.  There 
parameters  are  appropriate  to  the  pump  generated  by  a 
typical  laboratory  electron  beam  passing  through  an  argon 
plasma. 


b).  The  same  plasma  as  above,  but  with  a 
beam  which  is  now  faster  by  a factor  of  2.5.  The  pump  this 
beam  generates  is  here  a beam  mode  with  kQ  ■ 0.1. 


128 


to  a resonance  of  the  Stokes  wave  at  - k and  a low 
frequency  response  which  cannot  be  treated  by  a fluid  ap- 
proximation. Since  it  is  no  longer  true  that  u>«kve,  we 
have  had  to  make  the  approximation: 

(47)  ■ 

In  addition,  we  have  used  the  ordering: 

Linder  these  conditions,  the  solution  to  (6)  has  a maximum 
growth  rate: 

3W,S.‘S 
u ' 3 MK'k)  (S-* 

The  SM  instability  has  disappeared  because  of  the  heavy 


damping  of  the  anti-Stokes  mode. 

To  understand  why  Quon,  et  al.  do  not  observe  this 
OTS-like  structure,  we  must  study  the  narrow  beam  approx- 
imation, equation  (9),  on  which  our  application  of  (6)  to 
a beam-produced  spectrum  is  based.  Indeed,  it  appears  that 
this  reqion  of  growth  only  exists  in  the  narrow  beam  case. 
While  the  decay  region  easily  satisfies  the  narrow  beam 
condition,  for  the  OTS  instability,  we  have: 

)W-S_|  « 3kAk,  ( 

or 

uj  >>  00b  • 

Since  this  is  not  the  case  for  our  parameters,  the  broad 
beam  has  washed  out  the  OTS  instability,  so  degrading  the 
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growth  rate  that  Quon,  et  al.  do  not  observe  it. 

Two  changes  occur  in  this  electron  beam  experiment 
when  the  beam  velocity  is  increased,  as  in  the  second 
experiment  by  Quon  et  al.: 

1) .  kQ  becomes  smaller.  This  is  clear  from  the 
fact  that  the  fastest  growing  mode  of  the  beam  instabil- 
ity has  wavenumber: 

« un  vk/vb  • 

2) .  There  are  electrostatic  accelerations  which 
cause  the  beam  spread  to  decrease.  This  not  only  provides 
a relatively  cold  beam,  but  also  generates  a high  frequency 
longitudinal  wave  with  non-Langmuir  dispersion,  a beam 
mode  as  described  above.  For  Wong's  parameters,  based  on 
O'Neil  and  Malmberg's  results  for  a narrow  beam,  we  get: 

W,  ^ 1.05  7 , ko  -0.  I. 

We  have  been  able  to  treat  the  linear  stability  of  this 
wave  with  equation  (6),  by  merely  using  the  above  relation 
for  uie(k0)  instead  of  the  usual  Langmuir  dispersion.  No 
other  change  in  (6)  is  required,  since  its  derivation 
depends  only  on  the  pump  wave  being  high  frequency  and 
having  longitudinal  polarization.  The  numerical  solution 
of (6)  for  these  parameters  is  shown  in  Figure  eleven  b). 

This  figure  shows  two  parametric  processes,  a decay-like 
region  to  the  left,  and  an  OTS-like  region  near  kp  - k = kp. 
The  decay  region  is  not,  as  Quon,  et  al.  , state,  any  more 
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difficult  to  excite  than  in  the  case  of  the  slower  beam. 

In  fact,  since  the  matching  condition  can  be  satisfied  for 
somewhat  smaller  values  of  k0  - k,  the  Landau  dampinq  is 
smaller  and  the  growth  rates  for  this  case  slightly  larger 
than  for  the  previous  experiment.  Notice  that  perfect 
matching  occurs  when: 

I?)  = UJa(H)  ■ 

This  demands  that  k = 0.26  or  kQ  - k = -0.16.  The  growth 
rate  is  well  accounted  for  by  equation  (23).  Again,  (47) 
must  be  used  to  account  for  the  OTS-like  region  of  g-owth, 
although  here  we  have  the  ordering: 

CO  » fr(k*±k);  X.=0- 

so  that  the  solution  to  (6)  can  be  approximated: 

LO  Z il  /4  * 

since  here,  « 0.022.  Interestingly,  Quon,  et  al. 

report  neither  of  these  kinds  of  plasma  behavior.  Since 
these  are  the  results  of  a linear  relaxation  of  the  team- 
induced  spectrum,  we  have  concluded  that  the  effects  that 
Quon,et  al.  do  observe  are  due  to  non-parametric  processes. 

The  strongest  parametric  interactions  possible  are  those 
that  (6)  predicts.  That  they  are  not  observed  points 
towards  other  effects,  most  probably  trapping  effects, 
being  the  processes  that  Quon,  et  al.  observe  in  their 
second  experiment. 

. 

- » - 
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CONCLUSIONS 

We  have  developed  a theory  which  allows  us  to  treat 
the  problem  of  linear  relaxation  of  an  enhanced  Langmuir 
spectrum  without  approximation  in  k^j  and  taking  into  con- 
sideration multidimensional  effects.  Three  new  results  of 
note  come  out  of  our  analysis: 

1) .  In  a number  of  regimes,  it  is  possible  to  have 
the  Langmuir  pump  decay  into  daughter  waves  with  wavenumbers 
larger  than  the  pump  wavenumber.  Particularly,  at  equal 
temperatures,  the  SM  instability,  has  daughter  waves  at 

kQ  + k & kQ.  Thus,  these  waves  have  a phase  velocity 
slightly  greater  than  the  phase  velocity  of  the  beam.  These 
waves  may  be  able  to  accelerate  electrons,  and,  for  a 
beam-induced  pump,  draw  out  a tail  on  the  beam. 

2) .  In  every  regime,  there  is  significant  three- 
dimensional  structure,  which  will  lead  to  a spreading  of 
the  spectrum  in  angle.  The  full  problem  of  the  three- 
dimensional,  nonlinear  relaxation  of  the  spectrum  has  not 
been  investigated,  but  our  linear  results  show  that 

this  is  essential.  Especially,  when  the  OTS  instability 
is  off-axis,  large  growth  rates  can  occur  at  angles  exceed- 
ing 20°  to  kQ.  The  angular  spread  of  the  spectrum,  es- 
pecially in  the  study  of  relativistic  beams,  is  a critical 
parameter,  and  the  wave  effects  that  we  have  studied  here 
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affect  it  greatly. 

3).  For  a range  of  parameters  with  0e/Gi  y>  1 and 
kQ  in  a critical  range,  the  dominate  process  of  decay  is 
forward  scattering.  This  result,  coming  uniquely  from  a 
study  of  the  kQ  dependence  of  the  spectral  decay,  has  im- 
mediate bearing  on  especially  the  theories  relatinq  to 
Langmuir-to-electromagnetic  mode  conversion.  A whole  new 
set  of  kinematic  relations  is  possible  with  these  daughter 
waves;  the  k0  = 0 result  of  daughter  waves  with  equal  ana 
opposite  wavenumbers  to  the  pump  is  true  only  for  values 
of  kQ  large  compared  to  the  root  of  the  mass  ratio. 

To  be  sure,  our  results  are  only  qualitative  guides 
to  the  understanding  of  the  real  problem  of  Langmuir  tur- 
bulence. What  the  linear  theory  we  have  presented  here, 
offers  is  a physical  understanding  of  the  first  stages  of 
the  process  of  relaxation  of  a narrow,  Langmuir  spectrum  in 
a homogeneous,  field-free  plasma.  The  most  important  next 
step  in  this  study  is  to  extend  our  results  to  the  case 
with  a magnetic  field,  expecially  since  applications  of 
our  results  to  spectra  enhanced  by  electron  beams  frequent- 
ly occur  in  magnetized  plasmas.  And,  of  course,  the  full, 
nonlinear  problem  of  the  spectral  decay,  must  be  solved. 

Our  results  are  only  a beginning  to  that  study. 
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DERIVATION  OF  THE  SUSCEPTIBILITIES 


The  method  used  here  to  derive  the  susceptibilities 
is  due  to  unpublished  work  by  Martin  V.  Goldman.  We  Legin 
with  a fluid  equation  for  v as  a function  of  E;  all  other 


fields  and  densities  have  been  eliminated: 


(A.l) 


where : 


KK)  = roO  E(fc)  + J (£-  %r)-T(ic.)  v(Kt) • e(Kt) 

V(K,)-kt  [T(K)  + TCti'EQOl 

k l(K)  y(K,)  -HMUi 
+ S K T(K)  JdKmVCU  VCK,)  kyVCKj 


JKr 

d K123 


■d‘U*2$dU' 

d\,(Lj,U,ktd^(l\t4to i cy,  l,  l L\r/.  v 

(2n)"-  ® ~ -W,~  . 


The  strategy  is  now  to  calculate  the  various  terms  in  the 
expansion  of  the  velocity  in  terms  of  the  field,  use  that 
in  the  integral  for  the  current,  and  then  group  terms  by 
powers  of  E to  find  the  various  susceptibilities: 


(A. 2) 


i°W  = TiKbttK), 
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VW=  f2-^  frO:,)-£«)]{rOcym] 

Now : 

)»(K)  = -eJdKn  [nW(K,)V‘tl4)  * nw(K,)  vI,’(Ki)J 

(A.4)  * P'EW  ffe-TW-EdOj 

♦ TO'iVft^^  ItW  EOc.)]  . 


Comparing  this  with  the  definition: 

(A. 5)  j**(ti  = -iuifdKn  JCCK.K,, Ki)- £„(&)•%(&)  e, 
we  get: 

v/i/./  x ..  4mc’rw  f iklftoel.  \2^tj,(Ajt  2,  ■T(k)-Tt&)-2i  *>l 

r - 2nfufu,ak  1 eto  1 el 

- £ Tt*.)  etm  u^i  e,Tf«*)e*1 

(a  . 6 ) g e J 


_ w>xk  TftO  ■£,  et- ‘DK)<&  'W,  »»l  TlKi) dt- 1 (k.)  . £ 7 
e*  e*  J > 

and  since  two  of  these  frequencies  will  always  be  high,  we 
tave  for  any  instability: 


(A. 7) 


I bl.^V.  lo  l^L 


'Ut»Vt  Ki 


• - yutV/w1 

Ul%  U>  t ) 


Where  the  are  the  polarization  vectors  of  the  fields 
E(ki).  For  the  case  when  one  of  the  fields  is  a low 
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frequency  field,  only  one  term  contributes: 


(A, 8) 


c^nri  C ktj  7 r 7 . 

2m'ww,((*ttu))  Ck‘v*‘  ) l k,lk*k»l\ 


The  third-order  velocity  is  gotten  similarly  and 
consists  of  seven  integrals,  which  must  be  then  substi- 
tuted in  the  expression  for  the  third-order  current.  The 
integrals  which  are  dominant  for  the  case  of  two  high 
frequency  waves  are: 

f U Tft.).EC*a  g-lQ.T(K-K>)/fe-fet)  EQQ-fOO 
J ^ 1 4rh  IJ  U 0JX  (ui-cu,) 

t T(KOE(K,)  7k-^VTOKJ-(k->e.)  B(K 

Vi(w-CJi) 

+ rcK.yecK,)  fr-k .hik-k.)  ■»-*.)  f(ko-eck») 

Ld  X (_U>~Uli) 

+ T(K,1£(k,)  (k  fcl.RK-K.Xfe.-*.)  Elk) -EM 

ojiCuj-Wi) 

, TtlU-Hfc)  (k-fcOKii-KOft-fc)  KfeKflU 
4 m)'KK.)  (HaVTTK-lUft-fe)  £(fc) ■&>:■)  7 . 

4/(  (U/  -Ufi)  J 

Comparing  this  expression  with  the  definition  of  the  third- 
order  susceptibilities,  and  using  a travelling  wave  pump, 
we  find: 


(A. 9) 


4iri‘c,n* 


v)w  Y*1 


A ^ A 

e0*t  e.-e. 


and: 


X (K±K#,  K.,  *K#,  KH6) 


4rrie^ e.  e*  £•£. 


* (U0i  lo)  (w,  rto)  U4l  Ue-t 


(A. 11) 
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We  also  note  here  some  symmetry  properties  of 
these  susceptibilities,  which  can  be  derived  from  either 
the  integral  definitions  of  the  susceotibilities , or  from 
the  explicit  forms  we  found  for  them: 

(A*12)  //m*  (K.K,.Ki)  - 7-JIm  (*.  Kz.Ki) 

*0  r ( K,  Kh  Ki,  Kz ) 

(A. 13)  ) 

And,  using  the  these  relations: 


<A-14a)  tu.l-K-Kr&fimtoK.tyha  (h-K K), 

(a. 14b)  ■ ^i)  - Xdih/  K,  Kz)~ 


